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ABSTRACT 


Title of thesis: Topics in Astrophysical X-ray and Gamma Ray Spec- 

troscopy 

Roger Wo Bussard, Doctor of Philosophy, 1978 
Thesis directed by: Doctor Reuven Ramaty 

We have investigated a number of topics relating to astrophys- 
ical observations that have already been made or are currently planned 
of spectral features, mostly emission lines, in the X-ray and gamma 
ray regions of the electromagnetic spectrum. These topics include: 
the production of characteristic X-ray and gamma ray lines by non- 
thermal ions, spectral features induced by processes occurring in 
strong magnetic fields, and the positron annihilation line at 0.5 
MeV. Although astrophysical spectroscopy in this regime is quite 
young, a number of conclusions can be drawn from existing observa- 
tions and upper limits, and much can be learned from experiments 
currently under consideration. 

First, we have calculated the rate of X-ray production at 6.8 
keV by the 2p to Is transition in fast hydrogen- and helium-like 
iron ions, following both electron capture to excited levels and 
collisional excitation. We used a refinement of the OBK approxi- 
mation to obtain an improved charge exchange cross section. This 
and the corresponding ionization cross section were used to determine 
equilibrium charge fractions for iron ions as a function of their 



energy. The effective X-ray line production cross section was found 
to be sharply peaked in energy at about 10 MeV/amu. Since fast ions 
with similar energies can also excite nuclear levels, we have cal- 
culated the ratio of selected strong gamma ray line emissivities to 
the X-ray line emissivity. We use these calculations to set limits 
oh the gamma ray line intensity from the galactic center and the 
radio galaxy .Centaurus A, and we find that these limits are generally 
lower than fluxes reported in the literature. In addition, the 
presence of interstellar dust grains can give rise to a very narrow 
line profile for some gamma ray lines from nuclear excitations. 

In a grain, a nucleus with a long lived excited state may be stopped 
before emitting the gamma ray. We have compiled the relevant range 
and energy loss relations and have then shown the results of this 
effect for several examples. 

12 

The intense magnetic fields (of order 10 gauss) believed to - 
exist at the surfaces of X-ray pulsars could produce observable 
structure in their spectra. Such fields severely quantize the motion 
of electrons in the transverse plane, and line emission and absorp- 
tion can occur accompanying transitions between levels. We have 
confirmed a previous calculation of the rates for these processes, 
and we have calculated the cross section for electron-ion Coulomb 
collisions in strong fields. This cross section is significantly 
modified by the quantization due to the magnetic field. We have then 
calculated the line production rates for direct excitation by accreting 
material and by excitation from thermal, collisions . The results are 



compared with recent observations of Hercules X-l, and we find that 
the interpretation of the spectrum in terms of an emission line is 
less likely than that of an absorption feature. - Future observations 
with higher spectral resolution should resolve the question. 

Finally, we have studied the annihilation of galactic positrons 
in order to evaluate the probabilities of various channels of anni- 
hilation and to calculate the spectrum of the resulting radiation. 

The narrow width (FWHM less than 3.2 keV) of the 0,51 MeV line ob- 
served from the galactic center direction (leventhal et al. 1.978 a,b) 

implies that a large fraction of the positrons should annihilate in 

5 

a medium of temperature less than 10 K and ionization fraction greater 
than 0.05. HIX regions at the galactic center could be possible sites 
of annihilation. 
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I. INTRODUCTION 

As trophy sical information from our galaxy comes to us primarily in the 
form of electromagnetic radiation, from radio frequencies to high energy 
gamma rays, ‘In the attempt to understand the nature of various astrophysical 
objects and regions and their physical processes, much can be learned from the 
observations in particular of spectral lines and features (broadened lines, 
absorption edges, cutoffs, etc.) For example, red-shifted lines in the 
spectra of extragalactic objects lead us to believe that the universe 
Is expanding, siiice the amount of shifting is proportional to the distance 
to the object. Locally, lines from electronic transitions in atoms allow 
determination of relative elemental abundances, the spin flip line of atomic 
hydrogen at 21cm wavelength has been used to map neutral hydrogen in the 
galaxy, and the rotational transition lines of molecules give the' galactic 
distribution of dense molecular clouds.' Moreover, the simultaneous observation 
pf two or more lines usually provides information on the physical conditions 
in the source, such as the density and temperature. 

This work consists of a theoretical consideration of spectral features, 
mostly emission lines, from a very new astrophysical window, that of X-rays 
-and soft gamma rays (about 6 keV to 6 MeV), Some of the features considered 
have been observed and others are' possibly observable by currently planned 
experiments. Measurements in this spectral region must be carried out from 
balloons, rockets or satellites because of atmospheric absorption. The 
first line considered results from the' electronic transition from the 2p state 
to the Is state in highly ionized Iron at ~ 6.8 keV, the iron Ko; line. 

Elements with higher Z, and consequently higher Kor line energies, are far 
less abundant than iron, and lower energy K a line protons (from lower Z 
elements) will suffer much absorption, via photo ionization, propagating through 


1 



2 


the interstellar medium, The 2p state, ia populated foilwing charge exchange 

and radiative recombination to excited states and the subsequent cascade, or 

by excitation if the ion already has one or two- electrons-. Since these 

processes occur most readily when the relative velocities of the interacting 

particles is about one tenth of the speed of light, Xa emission can be produced 

‘ 7 

in a plasma with temperature ~ 10 K, or in regions with significant fluxes of 

ions with energies 10 MeV per nucleon. Tn the former case, the emitting 

iron atom is essentially at rest, and the line profile suffers little Doppler 

broadening. In the latter, the iron has the velocity, ~ ,1c, mentioned above, 

leading to a line width (FWHM) ~ 2.4 keV. Thus the broad 6.8 keV line is a 

tracer of low energy cosmic rays. Another type of line feature considered 

results from the excitation of nuclear levels in strong interactions by 

nuclei with the same energies. For example, a carbon 12 nucleus can be 

excited in a collision with an energetic proton, and it will subsequently 

deexcite by emitting a 4.44 MeV photon. Next, we consider the possibility 

that the strong magnetic fields associated with pulsars can give rise to line 

features analogous to the cyclotron lines from plasma in weaker fields.- -For 

12 

a field strength of 10 Gauss, for example, the fundamental cyclotron fre- 
quency corresponds to a photon energy of ~ 11 keV. Finally, we treat the 
production of the 0,511 MeV line from the annihilation of positrons with 
electrons for various galactic sites. All these spectral features have been 
Observed or could be observed by detector systems presently under consider- 
ation as for the Gamma Ray Observatory, for example. 

A search for the Ko; line of highly ionized iron was begun in the late 
1960*s and early 1970’s by rocket bom proportional counters (Holt et al. 

1968 and 1969) and Bragg crystal spectrometers (Kestenbaum et al. 1971 and 
Griffiths 1972). Although this search was not fruitful, the iron line 
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eventually was observed from some weaker sources with improved, more sensitive 

detectors: young supernova remnants (Serlemitsos et al. 1973, Brisken 1973, 

Becker et al. 1976, Pravdo et al. 1976, Davison et al. 1976), collapsed stars 

in binary systems (Serlemitsos et al. 1975, Pravdo 1976, Pravdo et al. 1977a), 

and clusters of galaxies (Mitchell et al. 1977, Serlemitsos et al, 1977).' 

Although it now appears that the line emission from young supernova remnants 

is consistent with that from a thermal plasma, low energy cosmic ray iron 

nuclei undergoing charge exchange with the ambient medium were originally 

proposed to explain this line (Serlemitsos et al, 1973). Nucleons with these 

« 

energies also produce gamma ray lines by exciting nuclear levels and in 
Chapter IX, we have calculated tire ratio of the gamma ray line emission to 
the iron Kor emission produced in this way. These calculations, together with 
X-ray observations are used to place limits on gamma-ray line production by 
fast ions in various sources . 

Astrophysical gamma ray spectroscopy is as young as that for X-rays . In 

two large solar flares of 4 and 7 August 1972, Chupp et al.(1973, 1975), 

using a Nal scintillation detector flown on OSO-7, reported observations of the 

0,5 MeV annihilation line, the 2.2 MeV line from neutron capture on hydrogen, 

the 4,4 MeV line of o and the 6.i MeV line of 0. Lingenfelter and 

hamaty (1967) had earlier made predictions of gamma ray line fluxes in solar 

flares, and interpretations of the observations may be found in Ramaty and 

Lingenfelter (1975), and Ramaty et al, (1975, 1977),' More recently, a balloon 

* flight bearing a high resolution, high purity germanium detector has observed 

a strong line at 510,7 HH3.5 keV from the direction of the galactic center with 

—3 —2 —1 ■ 

flux (1.21 +0.22)xl0 photons cm s (Leventhal, McCallum and Stang 1978a, 
b) , We consider the implications of this result in Chapter IV. In previous 
observations of the galactic center from two pioneering balloon flights in 
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November 1970 and November 1971, a group at Rice University, using a scin- 
tillation spectrometer, reported a marginal observation of a line at 476 +24 

-3 -2 -1 

keV with a flux of (1.8 +3.5) x 10 photons cm s (Johnson, Hamden, and 
Bayraes 1972, and Johnson and Haymes 1973). Despite the low significance of 
the feature due to background problems, many interpretations of this emission 

7 

were made, including; (i) cosmic ray excitation of interstellar Li to its 

\ t ? i 

478 keV level via \i(p,p' )^Li* (Fishman and Clayton 1972) and Sle(a,p)\i* 
(Kozl<?vsky and Ramaty 1974) , (ii) gravitationally redshifted positron annihil- 
ation radiation from neutron star surfaces (Ramaty, Bomer, and Cohen 1973 
and Guthrie and Tademaru 1973), and (lii) positronium annihilation from its 

triplet state via three photon decay (Leventhal 1973) . Then in April of 
» 

3,974, the same group, using a similar instrument but with smaller opening 
angle, failed to observed the above mentioned feature but saw instead a 

_4 —1 

feature at 530 +11 keV with a flux of (8 +3,2)xl0 photons cm s (Haymes 
et al. 1975). Aside from detector problems, time variability seems the most 
likely interpretation of the discrepancy. We note that in all three obser- 
vations, the significance of the features was ~ 3<x, 

The Rice University balloon flight in 1974 gave evidence for an emission 

line at ~ 4.5 MeV in addition to the 530 keV line from the direction of 

the galactic center (Haymes et al. 1975). We note that more recent observation 

from HEAO—1 has failed to confirm this feature, however. The ~ 4.5 MeV line 

was also seen from the radio galaxy Centaurus A, again at about the 3a level 

(Hall et al. 1976), A detailed treatment of the production of gamma ray lines 

by energetic particle reactions is given in Ramaty, Kozlovsky and Lingenfelter 

0.979); essentially all possibly observable lines are considered, and they 

12 

find that the line at ~4.4 MeV from the first excited level of C should 
he the most easily detected line produced, especially for the lower resolution 
Nal scintillation detectors. Observations of the gamma ray lines produced 
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by energetic particle reactions, together with the iron line discussed above, 

O' 

would incicate concentrations of low energy cosmic rays, and thus perhaps 
provide information on the sources of the cosmic particle radiation. 

Finally, we note that Jacobson et-al. (1978) have reported observations 

of a transient event in June of 1974 which included lines at 0.413, 1.79, 

2.22 and 5.95 MeV, seen by a Ge(Li) high resolution detector. An interesting 

interpretation of this event is given by Lingenf elter, Higdon and Ramaty 
. •>« 

(1978), in terms of gravitational redshifts. 

In Chapter II the production of certain characteristic X-ray and gamma 
ray lines by energetic but subrelativistic ions is considered. In the first 
section, we have calculated x-ray production at ~ 6.8 keV by the 2p to Is 
transition in fast hydrogen- and helium-like iron ions, following both 
electron capture to excited levels and collisional excitation. We used a 
■refinement of the OBK approximation to obtain an improved charge exchange • 
cross-section. This, and the corresponding ionization cross-section were 
used to determine equilibrium charge fractions for iron ions as functions 
of their energy. The effective X-ray line production cross-section was found 
to be sharply peaked in energy at about 8 to 12 MeV/amu. Since fast ions of 
similar energies can also excite' nuclear levels, we have calculated the 
ratios of the 4.4 MeV and 0.85 MeV (the most intense nuclear line of **^Fe) 
gamma ray line emissivities to the X-ray line emissivity. We used these 
-calculations together with X-ray observations of the galactic center from 
QSO-8 and HEA0-1 to set limits on the intensity of gamma ray line emission 
from the galactic center and the radio galaxy Centaurus A, and we find 
that these limits are generally lower than those reported. In the second 
section of Chapter II the effect of Interstellar dust grains on the width 
of certain gamma ray lines produced by fast ions is considered. Specifically, 
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for some relatively long-lived nuclear levels, the recoiling excited nucleus 
may be stopped in the grain before deexciting, which leads to a very narrow 
Tine profile. This problem was originally considered by Lingen'felter and 
Ramaty (1977); they have treated the 6.1 MeV line of ^0 in detail. One of 
the difficulties they encountered was the lack of consistent range relations 
through the relevant energy region (~ 1 to 100 keV/amu). In addition, since 

the range of the recoiling ion is of the same order as typical grain radii, 

* 

straggling (statistical fluctuations in range) is expected to be important. 
Therefore, for the calculations presented here, we have compiled a consistent 
set of range and energy loss relations, taken the effect of straggling into 
account, and considered lines from less abundant elements. Using a Monte 
Carlo simulation, we have evaluated the expected strengths of these very 
narrow lines as functions of the grain size distribution for various 
energetic ion spectra. Thus possibility of obtaining information from 
. these features with high resolution solid state gamma ray detectors is 
discussed. ■ 

12 

The expectation of intense magnetic fields Gs 10 gauss) at the 
surfaces of neutron stars gave rise to consideration of another kind of 
high energy spectral feature, since for such fields, the fundamental 
cyclotron frequency lies in the hard X-ray band. Classically, it is well 
known that a magnetic field will force electrons to move on. a helix centered 
on a field line called the guiding center. As they orbit, they produce 
gyrosynchrotron radiation, lose energy as a result and consequently spiral 
inward toward the guiding center. If, however, the field is very strong, 
the orbits of the electrons become quantized; the gyroradius of the electron 
orbit and its perpendicular momentum take on discrete values . Radiative 
transitions between levels can occur and this process results in a quanti- 
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zation of the gyro synchrotron radiation. Line emission occurs accompanying 
transitions from higher to lower transverse momentum levels, and upward 
transitions take place with absorption of resonant photons. 

* Several observed X-ray sources, in particular X-ray pulsars, are believed 
to be spinning neutron stars with the field anchored rigidly to- the surface. 
They are powerful sources and are thought to derive their energy from the 

accretion, of matter from a companion star. Hercules X-l is the best studied 

* 

of these objects; it shows significant spectral changes during a portion of 
its pulse phase (Pravdo et al. 1977b and 1978). Since in this picture the 
pulse phase would correspond to the angle between the line of sight and the 
magnetic field lines, these spectral changes are evidence for the effect of 
the magnetic field on photon production, absorption and scattering processes. 
Furthermore, there is evidence for an emission feature near 60 keV in the . 

Her X-l spectrum (Trumper et al. 1977 and 1978, Coe et al. 1977), which has 
been interpreted by these authors as a cyclotron line (see also Daugherty 
and Ventura 1977) . 

Although details of accretion flow onto rotating, highly magnetized 
neutron stars are unknown, if the rotation is fast enough, any plasma which 
falls to the surface is generally assumed to be channeled by the field into 
a polar cap with essentially zero pitch angle (Basko and Sunyaev 1975 and 
1976, Eisner and Lamb 1976). The infalling ions lose their free fall energy 
by heating the atmosphere at the polar cap. The magnetic field quantizes the 
transverse energies of the electrons, and transitions between levels give 
rise to line radiation. The levels can be excited by collisions with the 
infalling ions or by thermal collisions with ambient ions; however, the 
correct cross sections for Coulomb collisions in these circumstances was 
not known. So in Chapter III, we have calculated these by using a 
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relativistically correct quantum electrodynamic treatment; then we have used 

them to calculate the energy loss rate for fast zero pitch angle photons. 

It was found that this rate is- slower than that for- -a non-magneti-zed plasma, 

but it is not reduced as much as had been predicted (Basko and Sunyaev 1975). 

Including the effect of elastic proton— proton scattering, we have calculated 

the multiplicity of cyclotron line photons produced by an ion falling onto 

a highly magnetized, hot plasma, and the fraction of its energy that goes 
* 

directly into line photons. We have also calculated the rate of excitation 
of the cyclotron levels by thermal electron-ion collisions. These results 
are presented in Chapter III, where we also discuss their relevance to 
observations of Her X-l. 

Strong magnetic fields will also affect the X-ray spectra of the plasma 
directly, by both modifying the Compton scattering cross section (Canuto 
et al. 1971, and Gnedin and Sunyaev 1973) and giving rise to resonant 
f cyclotron absorption, i.e„, absorption of a photon accompanied by an exci- 
tation of the electrons transverse energy (Daugherty and Ventura 1978). 

An attempt has been made by Boldt et al, (1976) to explain the high energy 
0- 20 keV) steepening observed for Her X-l (Pravdo 1976) by Compton scattering 
off thermal electrons using the cross section of Canuto et al. (1971). 

However, in Chapter III, we have shown that cyclotron absorption could be 
important in forming this feature; specifically, we have found that the 
infalling ions can produce a .flux of relatively high, energy (~ few hundred 
keV) . electrons streaming down the field lines toward the stellar surface. 

The velocities of these electrons are large enough that they can absorb 
photons with energies below the cyclotron energy because of the Doppler shift 
of the photon's frequency. The production rates of the highen energy electrons 
and their effect on the x-ray spectra of the pulsars is discussed in 
Chapter III. 
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In Chapter IV, we have considered the formation of the 0.511 MeV line 
by positron annihilation In the galaxy. As was mentioned above, a narrow, 
intense line has been seen from the galactic center direction at 0.5107 
+0.0005 MeV (Leventhal et al. 1978a, b) with a full width at half maximum 
less than the detector resolution, about 3.2 keV'„ Since previous estimates 
of the line width have been larger (e.g., Stecker 1969), in Chapter IV, 
we have studied the annihilation of galactic positrons in order to evaluate 
the probabilities of various channels of annihilation and to calculate the 
spectrum of the resulting radiation. These channels include positronium 
formation by charge exchange with neutral atoms and by radiative recombina- 
tion with. free electrons and direct annihilation with either bound or free . 
electrons. In a nearly neutral gas, charge exchange will occur before most 
of the positrons can thermalize; otherwise they thermalize before annihilatiiig. 
We have calculated the fraction which annihilate in flight (before thermal i- 
zation) as a function of the ionization fraction of the medium by a Monte 
Carlo simulation, and we have evaluated the dependence on gas temperature 
of the annihilation rates for the various processes mentioned above for 
those which thermalize. Then, we have calculated the annihilation radiation 
spectrum for all these cases. 

We find that the narrow width of the 0.511 MeV line observed from the 
Galactic Center implies that a large fraction of positrons should annihilate 
in a medium of temperature less than 10"*K and ionization fraction greater 
than 0.05. Since most of the matter in the galactic center region is 
expected to be in neutral molecular clouds, it is difficult to reconcile 
the narrow line observation with diffuse positron production unless they 
are unable to penetrate the clouds. Another possibility is that the 
positron source is localized very near to the galactic center itself, in 
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which case the HII regions at the nucleus are possible sites ox ami j.uxj.ation . 

If future observations- confirm the small angular extent of the source of 
this. radiation- a massive collapsed object at the galactic center is -suggested. - 
Chapter V is a summary of the preceding topics, with a discussion- of 
possible future investigations, both theoretical and experimental. Searches 
for the characteristic iron X-ray line and nuclear gamma ray lines will 
provide constraints on cosmic ray origin- theories. High resolution spectros- 
copy 'can be useful in several ways. Nuclear lines from interstellar dust 
could be distinguished, depending on the size distribution of grains. 

Physical processes in the intense magnetic fields of binary X-ray pulsars 
can be determined by such detectors made sensitive to as low as ~ 30 keV. 

And high resolution spectroscopy at 0.5 MeV can provide information on sites 
of positron annihilation. 



II. X-RAY AND GAMMA SAY LINE PRODUCTION BY NONTRERMAL IONS 
a) Introduction 

A flux of fast subrelativistic ions in a medium can produce line 

emission in both the X-ray and the y-ray regions,, This chapter is concerned 

with the production of the lines of the 2p to' Is transitions in hydrogen- 

and helium-like iron ions (at ~ 6.8 keV) , and with the comparison of these 

lines with the strong y-ray line at 4.44 MeV resulting from the nuclear 
' 12 

deexcitation of C. The 2p to Is transitions in more abundant species 
such as C, N ,nd 0 give rise to less energetic X-ray photons, which subse- 
quently suffer interstellar absorption. 

The X-ray calculations presented in this chapter are of the same 
nature as those of Pravdo and Boldt (1975) who have treated X-ray line 
production from energetic oxygen nuclei in the interstellar medium. Watson 
(1976) has calculated the multiplicities of X-rays produced by fast iron and 
oxygen nuclei interacting with ambient hydrogen, but as we shall show, inter- 
stellar oxygen and iron, and not hydrogen, are the main contributors to elec- 
tron capture by fast iron ions. By using improved cross sections for charge 
exchange and excitation, we have carried out a detailed calculation of 
X-ray line production from fast ions moving through an ambient medium with 
solar system abundances. 

Astrophysical y-ray line production has been investigated by 
Ramaty, Kozlovsky and Lingenfelter (1975) and Ramaty, Kozlovsky and Suri 
(1977) for solar flares, and by Meneguzzi and Reeves (1975), Lingenfelter 
and Ramaty (1977), and Ramaty, Kozlovsky and Lingenfelter (1979) for the 
interstellar medium. By combining these results with the calculated X-ray 
line emission, we evaluate ratios of the X-ray and y-ray line intensities. 
These can be used to set limits on the expected y-ray line emission from 
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/arious source regions for which X-ray observations are available- As was 
nentioned in Chapter II, iron line emission has been observed from young 
supernova remnants and clusters of galaxies (Virgo , Coma, and Perseus), but 
this emission has a narrow width, consistent with production in a thermal 
plasma with a nearly solar abundance of iron. A broadened iron feature has 
been seen from. collapsed stars in binary systems, but the broadening is most 
likely due to rotational bulk motions, since charge exchange by the accreting 
matter is not efficient enough to provide the observed line luminosity. So, 
only upper limits are available at the present time on the intensity of the 
~ 6.8 keV lione from non thermal particles, 
b) X-ray Production 

1) Charge Exchange, Excitation and Ionization Cross Sections 

We consider energetic (about 1 to 300 MeV/amu) Fe ions 
moving through an ambient medium with solar abundances (Cameron 1973, except 
for the abundance of He for which we use the results of Brown and Lockman 
1976), The 2p to Is transitions in such ions having one or two bound 
electrons give rise to Ka X-rays at 6,9 and 6.7 keV, respectively. Because 
of Doppler broadening, however, these 2 lines marge into a broad feature 
centered at about 6.8 keV. For bare nuclei, the 2p state can be populated 
by charge exchange with ambient atoms; for ions having one electron, the 2p ' 
state can.be populated by charge exchange or excitation; and for ions having 
2 electrons, the 2p to Is transition can only occur following excitation. 

The method for calculating the charge exchange cross 
sections for bare Fe nuclei, a refinement of the Oppenheimer, Brinkman, 
and Kramers (OBK) approximation; is described in Bussard, Ramaty and Omidvar 
(1978). The OBK method is a Born approximation which treats target electrons 
as distinguishable (Oppenheimer 1928, Brinkman and Kramers 1930). Nikolaev 
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(1967) has developed an empirical factor which brings the OBK cross sections 
into agreement with existing data; we have applied this correction to the case 
of iron projectiles. The results at five energies, tabulated by the species 
of target atom, are shown in Table IX— 1. The three rows for each target 
element give the capture cross sections to any state which subsequently - 
deexcites to the Is, 2p and 2s levels, respectively. 

The charge exchange cross sections which lead to electrons 
in the 2p level are shown in Figure II-l as functions of the temperature of 
the ambient medium for the 5 energies given in Table 1. At low temperatures 
C*' 5x1 O^K) these cross sections are the same as those given in Table 1. The 
cross sections decrease with increasing temperature because the target atoms 
lose some of their orbital electrons. We have estimated the temperature 
dependence shown by folding together the contribution of the partial cross 
sections for charge exchange from each electronic orbital and the ionization 
fractions at a given temperature. For these partial cross sections we used 
the calculations described in Bussard, Ramaty and Omidvar (1978) , while the 
ionization fractions are from J. C. Raymond (private communication 1977). 

The results indicate that the cross sections are fairly constant up to a 
gas temperature of a few hundred thousand degrees, hut fall off rapidly with 
increasing temperature. It should be noted that the temperature dependences 
shown in Figure II-l might require some modification since the partial cross 
sections were obtained for neutral targets. In subsequent calculations, 

* however, we only use the cross section at low temperatures where the targets 
which are the major contributors to charge exchange are mostly neutral. This 
cross section is shown as a function of energy by. the dotted line in Figure 


II-2. 



The cross section for excitation of a hydrogen-like ion 


from the Is to the 2p configuration is given by: 

ivi *«rr <*■* 

(Bates 1962). Here Z is the atomic number of the gas atom causing excita- 

m 

tion, A is the Bohr radius (.529 a)» and Z is the charge, in units of the 
o 

electron charge, of the ion being excited. Also, m g and m^ are the electron ' 
and prpton masses, respectively, ,E is the kinetic energy per nucleon of the 
projectile, and 1^ is the ionization potential o£ hydrogen (13.6 eV). The 


lower limit on the integral, t . , is given by 

mm 
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where AE is the difference in energy between the two levels, and M is the 


reduced mass for the encounter. 


The ionization cross section of iron ions in collisions 
with ambient neutral hydrogen were taken from Watson (1976), but since he 
considered a gas consisting only of hydrogen, we have scaled his cross 
sections to account for the effects of helium. There are two contributions 
to these cross sections: one in which the target atom is not itself excited* 

called the elastic part, and an inelastic part, which includes the sum over 
all possible excitations of the target (e.g. Rule 1977). For iron projec- 
tiles on hydrogen and helium, each part contributes almost equally, but the 
elastic part scales as the target charge squared, while the inelastic part 
varies linearly with the target charge. Also we compared the calculations 
of Rule (1977) with the Watson (1976) cross section for ionization of the 
hydrogen-like iron and found good agreement. 

Since the energy required to ionize the iron projectiles 
is so large compared to the target ionization potentials, the ionizing 
collisions take place at very small impact parameters. Thus, electronic 
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screening of the target nuclei is not very important, and these same results 
should apply to an ionized medium. In fact. Rule (1977) calculated the 
ionization cross sections for both cases, and found differences of < 10% for 
the energies under consideration. 

2) Effective X-ray Production Cross Sections 

The effective X— ray production cross sections were 
obtained by folding the cross sections for capture and excitation with the 
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appropriate charge fractions: 
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where cr ^ (2p) is the effective cross section for X-ray production following 

the capture of an electron which cascades through the 2p to Is levels, 

<j c 1 (2p or 2s) is the cross section for the capture of an electron into 

the 2p or 2s level, either directly or following a cascade, by an ion which 

already has n electrons, a is the X-ray production, cross section following 

©xc 

excitation, and <7 is the cross section for excitation of an ion with 

exc 

n electrons. The cross sections cr c ^ n ^ are obtained by multiplying the cross 
sections given in Table 1 by (Z^^/26)'’, where is the bare charge corrected 

for electronic screening. The screening effect was taken from Bums (1964) 
for all the cross sections considered. The 2s term represents the magnetic 
dipole transition, which occurs more rapidly than two photon emission only 
in the triplet state of the heliumt-llke ion. On the average, this state is 
populated 3 times more often than the singlet state, giving rise to the 
factor 3/4. The f are the charge fractions or the fractions of all ions 


with n electrons 





The charge fractions satisfy the equation: 

+ 3'E O'r. I)'’ (t, T"' f } nti ^lN,) , i^ U 0 : (jr-4) 

j where the subscript I refers to ionization, is. the hydrogen density,, y 

dE 

the projectile velocity, and is the time rate of energy loss. In a 
steady state situation and if the energy losses are slow compared to charge 


exchange and ionization .collisions, equation II-4 reduces to 
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For this equation to be valid, the time between collisions represented by 
terms on the right side of equation II-4 should be much shorter than the 
energy loss time of the ions. Since the. charge exchange collisions manifest 
themselves in X-ray production, we expect that this condition will be satis- 


fied when, the X-ray multiplicity, defined by 


. n» r, p< 

(J- Y g 

is much larger than unity. . Here p is the gas mass density. 




a ' - a (2p) + a + a (2s) is the total 6.8 keV X-ray production cross- 

a C vXC o 

section from equation II-3, dE/dX is the energy -per amu lost per g cm of 
- matter traversed, and E is some energy at which the ions are essentially 
bare. In evaluating equation II-6 we assume that the Fe ion is in charge 

• equilibrium at all energies. The energy loss rate dE/dX in a neutral medium 

- ->. • 

was taken from the tabulations of Barkas and Berger (1964) and Northcliffe 
and Schilling (1970), for a gas composed of H and He with 8.5% He by 
. number. The effective charge for the Fe ions, was taken from Northcliffe 
and Schilling (1970). The losses of fast ions in an ionized medium were 


calculated by using 
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.Ginsburg and - Syrovatskii (1964), where e and m are the electron charge and 
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mass, Ze the ion charge, v is the velocity of the ion, and n is the electron, 
density, taken equal to the hydrogen density. . 

Table 11-2 shows the multiplicity of ~ 6.8 keV X-rays as 
evaluated from equation 11-6. for several values of E. For a neutral medium* 
the multiplicity is high enough that equation Ii-5 should be a good approxi- 
mation. However; in an ionized medium the multiplicity is low, and hence 

' I 

equation 11-5 may not be a good approximation. An exception is the case of 
continuous acceleration where another term is added to equation 11-4 to 
cancel the effects of the energy losses. But even if such acceleration is 
not present, we would still expect about one X-ray photon per each energetic 
iron ion. Since the multiplicity for the ionized case in Table 2 is already 
close to unity, we do not expect that the use of equation II-5 will lead to 
any substantial error, even in this case. 

Figure II- 3 shows the charge fractions as functions of 
the projectile kinetic energy per nucleon, calculated by using equation 11-5. 
It shows that on the average, the iron nuclei are bare above ~ 40 MeV/amu, 
thqt they have picked up one electron by ~ 11 MeV/amu, that they have two 
electrons by ~ 10 MeV/amu. Below about 6 MeV/amu, the average iron ion has 
3 or more electrons. When this happens, the outer electron usually takes 
the energy from collisions, so the Kor transition occurs rarely from ions 
with 3, or more electrons. 

The charge fractions shown in Figure II-3 were then used 
to weight the cross-sections discussed in section 11a, according to equations 
II-3. The results are shown in Figure II-2 as the solid lines. Both curves 
show the same energy dependence, and excitation is seen to be the dominant 
mechanism for X-ray line production. Both cross-sections are fairly sharply 
peaked in energy at around -10 MeV/amu. Because of the sharpness of the peak. 
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we estimate the width of the broad 6.8 keV feature to be =*• 2.4 keV. Also, 

.observation of this line will give the amount of iron with energy from approx-' 

-44 

imately 5 to 20 MeV/amu in the source. This number is about 10 of 

_3 

iron per photon per sec produced, for a hydrogen density of 1 cm , and 
varies little with the assumed spectrum of ions (+ ~ 50% over the spectra 
considered in this chapter) . 
c) Results 

We have integrated the X-ray production cross-sections shown 

in Figure II-2 over a variety of spectra of fast ions. The abundances -of 

these ions were extrapolated from measurements at > 300 MeV/amu (e.g. .Meyer, 

Ramaty, and Webber 1974), H: He: C: N: 0: Fe = 1: .1: 5.6x10 3 : 2.8X10 
-3 -4- 

5.4x10 : 5x10 , and taken to be independent of energy. The spectra used 
for all species are power laws in kinetic energy per amu with a low energy 
cutoff, E c : 

'K CE') 

where n(E) is the number density of ions per interval of kinetic energy/anru. 

We have ‘also set n(E) equal to zero for E>300 MeV. As mentioned above the 
gas abundances used were taken from Cameron (1973), with the exception of 
He (Brown and Lockman 1976). They are given in Table XX-1. 

The solid curves in Figure II-4 show the 6.8 keV line emissivity 

_3 

per H atom for an energy content of 1 eV cm in fast ions from 1 to 300 
Hev^amu. If can be seen that the spectra which maximize "X— ray production 
have breaks at around 10 MeV/amu and fall off rapidly above that energy. 

This is a result of the sharp maximum in the X-ray production cross section 
at around 10 MeV/amu, shown in Figure II-2. Because of evidence that a 
galactic gradient of heavy abundances exists, we have also repeated these 
calculations for the same fast ion fluxes in a medium with these abundances 
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increased by about a factor of 5, as Ramaty, Kozlovsky and Lingenfelter -(1979) 
have assumed for the 5kpc ambient medium based on the radial gradients found 
by Peimbert et al. (1978) for N and 0. These abundances are given in 
parenthesis in table II-l. The results are shown by the dotted lines in 
Figure II-4. The overall X-ray production is not affected much; however, due 
to the enhanced possibility for charge exchange, the production is shifted 
slightly to emphasize higher energy particles. 

Figure II- 5 shows the ratio of 6.8 keV photons produced to the- 
energy deposited in the gas by the fast ions, assuming local abundances. 

Here, the X-ray line emission is maximized with respect to energy deposited 
when the spectrum has a break at about 15 MeV/anu. The curves are shown for 
both a neutral and an ionized gas. They have the same shapes for both types 
of medium, but since energy is lost faster in a plasma owing to collective 
long range interactions (emission of plasma waves) the curves are displaced 
downward in an ionized gas n 

In addition to the X-ray lines of iron, ions in the energy range 

under consideration will cause nuclear reactions, with subsequent y-ray 

line emission. Meneguzzi and Reeves (1975) and Ramaty, Kozlovsky and 

Lingenfelter (1975) have shown that in these reactions the strongest line 

12 

is at 4.44 MeV from the first excited nuclear level of C. This line has 
two components. One is narrow, resulting from nuclei in the ambient gas 
excited by fast H and He, and the other is broad, from the excitation of the 
fast nuclei by ambient H and He. Since we are assuming different relative 
abundances of heavy nuclei in the gas and the fast particles, these two 
components have equal intensities, and hence are treated separately in the 
present work. For the calculation of both the narrow and broad 4.44 MeV 
line intensities, we use the cross section compiled by Ramaty, Kozlovsky 
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and Suri (1977). In additions we also Include the broad component of the 
0.85 MeV line of. ^Fe, since this line intensity scales directly with the 
amount of fast iron which produce the broad ~ 6.8 keV X-ray line.. . 

In Figure II- 6, we show ratios of y-ray to X-ray line emissivi- 
ties, qy/q(6.8), calculated for the equilibrium case, for the two components 
of the carbon-line and the broad component of the iron y-ray line. In addi- 
tion to the effects of direct excitation of carbon by hydrogen and helium, 

we have also included the contribution of spallation of nitrogen and oxygen 
\ 

to the excited carbon nucleus. The X-ray production scales directly with 
the relative- abundance of iron in the fast ions, but varies in a more compli- 
cated way with the relative abundances of heavy nuclei in the gas, primarily 
iron and oxygen. On the other hand, the broad- component of y-radlation 
scales with the relative abundance of fast C, N and 0 ions, and the narrow 
component scales with the gas abundances of these nuclei. Thus, the ratio of 
y-ray to X-ray yields depends somewhat on the relative abundances of medium 
and heavy nuclei. For example, by increasing the abundances of C and heavier 
nuclei in both the gas and fast particles by a factor of 10, the 4.44-MeV 
y-ray intensity also Increases by a factor of 10, but the ~ 6,8 keV X-rays 
.increase by about a factor of 20 at' most ^ and hence q^/q(6.8) remains constant- 
or decreases somewhat, up to a factor of 2. 

Figure II-6 shows that there is less than an order of magnitude 
of variation of q^/q(6.8) with the various spectra assumed for the fast ions, 

except for those which flatten at ^ 50 MeV/amu. The flatter spectra favor 

✓ 

y-ray emission because the cross sections for nuclear excitation drop off 
much more slowly at high energies than does the X-ray production cross section. 
This trend is also reflected in the fact that lower spectral indices favor 
y-ray emission. Since we have assumed a high energy cutoff at 300 MeV/amu, 



21 


we do not take into account y~ ra y emission from tt° decay produced by particles 
at higher energies. 

d) Applications to the Galactic Center, Cen A, and Cas A 

We now apply the results of the previous section to examine the 

consistency of the reported 7 - ray data from the galactic center region and 

the radio galaxy Centaurus A (Haymes et al. 1975, Hall et al, 1976) with 

upper limits on the nonthermal 6.8 keV line from these sources (Kellogg et al. 

1971, R. Mushotzky, private communication 1977). We also set limits on the 

flux of energetic particles in the supernova remnant Cassiopeia A by using 

upper limits' on the broad X-ray line emission from this source (P. J. 

Serlemitsos, private communication). 

Haymes et al, (1975) have claimed to see a feature at around 

4.6 MeV from the direction of the galactic center which they attributed 

12 

to the 4.44 MeV line of C. The opening angle of the y-ray detector used 

•was 13° (FWHM). From the observations of Kellogg et al. (1971) of the X-ray 
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sources GCX at the galactic center, we take the measured flux of 1.1x10 
- 2-1 

photons cm s in a 2.4 keV band centered at 6.8 keV as an upper limit on 
the intensity of the broad iron Ko; line from GCX. Using this upper limit, 
we have set upper limits on both the broad and narrow y~ray lines produced 
by fast ions in GCX. Under the most favorable conditions for Y - ray produc- 
tion considered (s = 1.2, E - 100 MeV/amu, see Figure 6 )- the upper limits 

-5 -4 -2 -1 

on the narrow and broad components are 3.2x10 and 3.3x10 photons cm s , 

-4 

respectively. Haymes et al. (1975) give an intensity of (9.5 +2.7)xl0 
photons cm s with a width of several hundred keV, so they could have 
been observing the broad component. In this case, our upper limit is a 
factor of three below their observation. Since the opening angle of the 
y-ray detector is larger than the angle subtended by GCX, Worrall et al. 



(1978) have examined data from the Goddard Cosmic X-Ray Spectroscopy Experi- 
ment on board 0S0-8 to obtain the X-ray flux from an area of 20°x20° centered 

on GCX„ They find the upper limit to a broad 6.8 keV feature is 0.13 photons 

-2 -1 -3 
cm s , which implies after scaling to 13°xl3°, upper limits of 8x10 

8x10 photons cm “s to the broad and narrow y-ray components, respec- 
tively. It should be noted the region is source confused, and efforts are 
presently under way to eliminate problems from discrete sources. Thus the 
observation of Haymes et al. (1975), provided they indeed represent the 
broad y-ray component and their source is larger than GCX, is not inconsistent 
with the upper limits set by the X-ray data. However, as we shall now discuss 
the energy density in energetic particles required to account for the y” ra y 
data under the constraints of the X-ray upper limits is much larger than 
the limits set by the distribution of gas perpendicular to the plane in the 
galactic center region. 

The energy density in fast ions required to produce a ~ 6.8 keV 

X-ray flux corresponding to the. gamma ray observation for s = 1.2 and = 

•100 MeV/amu, of 4.8x10 ^ photons cm ^ is, from Figure II-4, 4 x10 1 VMjj eV 

cm , where is the mass o.f hydrogen in the source, in solar masses. 

Scoville, Solomon, and Jefferts (1974) and Mezger (1974) have estimated that 

7 8 

the amount of gas in the nuclear region is some 10 ' to 10 M^, and so, the 

-3 

•above •y-ray flux requires several keV cm in fast ions. This energy density 
.is 2 orders of magnitude larger than that required to maintain the z distri- 
bution of gas at the galactic center (Sanders and Wrixon 1973). We could 
lower this energy density in a manner similar to that suggested by Ramaty, 
Kozlovsky and Lingenfelter (1979). They proposed that the abundances of C 
and heavier nuclei relative to H in both the ambient gas and energetic 
particles could be larger than solar and cosmic ray abundances. Also, we 
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can assume spectral parameters for the energetic particles which maximize 
the y— ray line emission with respect to energy content (s >3, 15 MeV/amu) 

• However, both these assumptions lower q^/q(6.8). As can be seen from 
Figure 11-6, for s = 3 and E^ = 15 MeV/amu, q^/q(6.8) is lower by about an 
order of magnitude than for s = 1.2 and E c = 100. MeV/amu. Furthermore, since 
q(6.8) depends on the relative abundances of the heavy nuclei in both the gas 
and the energetic particles, while q^ depends on these relative abundances 
in only one component, an increase of heavy nuclei abundances decreases 
q^/q(6.8). Therefore, we conclude that it is difficult to lower substan- 
tially the energy density required to produce the y-ray flux observed by 
Haymes et al. (1975) without leading to conflict with the X-ray data. The 
only exceptions could arise if the source region is so hot that the charge 
exchange and X-ray production are suppressed, or that it is opaque to the 
' X-rays but not to the gamma rays „ 

We note, from Figure 11-5, that to produce a 6.8 keV X-ray flux 
-2 -2 -1 

of 4.8x10 photons cm sec , the energy deposited in a neutral gas by the 

energetic particles with spectrum s = 1.2 and E £ = 100 MeV is about 
42 

4x10 erg/sec. This energy release is large, but still of the same order 
as the observed far infrared luminosity of the galactic center (Hoffman, 
Frederick, and Emergy, 1971). The energy deposited in an ionized gas is 
larger by a factor of 5 than that in a neutral gas. 

I If, nonetheless, we use the spectral- parameters s = 3 and 

- E c = 15 MeV/amu, we find that the observed upper limit on the broad 6.8 keV 

•line from a 13°xl3° region implies, from Figure II-6, upper limits of 
—5 —4 —2 —1 

: 1.4x10 and 2.8x10 photons cm s on the narrow and broad 4.44-MeV 

tline intensities. For these fluxes the energy density in fast ions is, from 

q 

(Figure II-4, only 7.4x10 /M„, a factor of 50 less than that required to 



LCcount for the observations of Haymes et al. (1975) under the other 
tssumed spectrum. 

One interesting possibility for explaining the incompatibility 

the gamma ray observation with the X-ray upper limits and gas distribution 

is time variability. In September 1977, observations of the galactic center 

from the HEAD-1 A4 experiment, a scintillation detector with opening angle 

-4 -2 -1 

20°, have provided an upper limit of 5.7x10 photons cm s in a 700 keV 

band centered on 4.4 MeV (J. Matteson, private communication, 1978). If we 

accept both measurements, then the galactic center source must vary on a 

time scale of a few years. Since this would imply a fairly localized source, 

we are tempted to suggest accretion onto a massive black hole as the source, 

as has been suggested to explain the time variability of Centaurus A (Fabian 

et al. 1976). Lingenfelter, Higdon and Ramaty (1978) have treated the y-ray 

line production in the accreting gas. Ions falling into the disk have little 

ox no time to cool, and predicted temperatures are favorable for nuclear 

interactions. This interpretation, however, has the disadvantage that the 

temperatures are high enough to suppress X-ray production by charge exchange. 

In addition to the galactic center, y-radiation was also detected 

from Centaurus A (Hall et al., 1976). They claimed a feature at 4.5 MeV, 

most likely less broad than ~ 150 keV, at a confidence level of 3.3cr, with 

—4 —2 —1 

a line flux of (9,9 +3.0)xl0 photons cm s . We examined the X-ray 

spectrum of Cen A around 6.8 keV (R. Mushotsky, private communication), and 

-3 -2 -1 

placed an upper .limit of 8.1x10 photons cm s on the broad Fe K& flux. 

Again, using the most favorable spectrum for y-ray production, we obtained 

-5 -4 -o -1 

upper limits of 2.5x10 and 2.6x10 photons cm “s to the narrow and 
broad 4.44 MeV components, respectively. Thus, the upper limit on the narrow 
component is a factor of 40 below the observation. Because of the quoted width 
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of the ~ 4 .5 MeV feature from Cen A, it is unlikely that this feature is the 

t 

broad 4,44 MeV line. 

Since supernovae or supernova remnants are candidates for sites 

of particle acceleration, we examined the spectrum of Cas A, a young supernova 

remnant. Because of a narrow thermal line at 6.7 keV (Pravdo et al, 1976), 

we used for an upper limit on flux in a broad line the flux above continuum 

at 4,85 keV (P, J. Serlemitsos, private communication 1977). This was fit to 

a gaussian centered at 6.8 keV and a full width at half maximum of 2.4 keV 

-3 -2 -1 

and gave an upper limit of 7x10 photons cm s in the broad line. At an 

assumed distance of 3 kpc, we obtain for the mass of Pe in the energy range 

5 to 20 MeV/amu, as discussed in section b2, <0,l/n u K , where n is the 

atomic hydrogen density at the source. From Figure II-6, we can set upper 

12 

limits on the broad and narrow components of the 4.44 MeV C line of 

-4 -5 -2 -1 

2.8x10 and 2.7x10 photons cm s , respectively. Presently planned 
y^ray. experiments could detect such fluxes, and thus we expect to obtain 
even more meaningful limits on the fluxes of energetic ions in supernova 
remnants from these experiments. 

e) The Effect of Interstellar Dust on Gamma Ray Line Shapes 

There is mounting evidence that significant fractions of the 
interstellar heavy nuclei are condensed out of the gas phase in dust grains 
with linear dimensions distributed around a few tenths of a micron. When 
these nuclei are excited by fast protons or alpha particles, they attain 
recoil velocities from around one keV per amu for ~*^Fe to several tens of 
keV per amu for CNO nuclei. In the gas phase, the nuclei are not decelerated; 
the recoil velocity is the velocity at which the gamma ray is emitted and 
consequently is responsible for Doppler broadening of the line. However, 
in grains, nuclei excited to states with long enough lifetimes have a chance 



to "be significantly slowed or even stopped before emitting the gamma ray, 
which leads to a much narrower line profile than produced in the gas, 

- To investigate, the importance of ..this effect quantitatively s we 

J 12- 16. 20 T 24„ 

have compiled the energy loss and range relations for C, O, Ne, ■ Mg, 

pQ Cd‘ ^*1 

Si, and °7e with energies from 10 to 1 MeWamu slowing down in water. 

t i. ^ 

For a recoil energy E greater than 12.5 keV/amu, these quantities are taken 

X *■ *' 

* % 

directly from the tables of Northcllffe and Schilling (1970). At lower 

energies, we have used the tabulated values of Winterbon (1975), using the 

scaling laws given there to apply his results to our nuclei and target 

material. Where overlaps occurred in the results from the two sources, 

differences were found to be less than about 10%. Another effect we had to 

consider was that of range straggling, or statistical fluctuations about 

the mean range. This quantity was obtained from the tables of Brice (1975), 

and the ratio of the variation in range to the mean range was found to be 

essentially the same for all nuclei considered. The rate of energy loss, 

dE/dX, the mean rectilinear stopping range, <X(E^)>, and the square of' the - 

ratio of the fluctuation in X due to straggling to X, (AX/X) are shown in 

Figure 11-7 for ^0, ^Mg and ~^Fe nuclei slowing down in water assuming 

that the grains are predominantly ice. These quantities, however, are 

not strongly dependent on the assumed grain composition. Fox a density 
—3 

of 1.5g cm , nuclei with recoil velocities as mentioned above slow down 

-13 -4 

in less than about 7x10 sec over distances less than ~ 10 -cm* By 

comparing with- the mean lives listed in Table II-3, (from Ramaty, Kozlovsky, 

and Lingenfelter 1978) and the expected sizes of interstellar grains, it is 

evident (Lingenfelter and Ramaty 1977) that many, lines could have significant 


very narrow components. 
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Ramaty, Kozlovsky and Lingenfelter (1979) have performed an 
extensive, detailed calculation of the gamma ray line shapes resulting from 
nuclear excitation hy energetic particles, using a Monte Carlo simulation 
which takes into account the reaction kinematics, excitation of the fast 
ions as well as the ambient gas nuclei, and excitation following spallation 
reactions, for various ion spectra parametrized hy power law indices and 
break energies as in equation II-8. In this calculation, a large number of 
interactions are considered; the interacting particles and the energy of the 
nonthermal ion are chosen according to the cross sections for the inter- 
actions, the abundances of the particles and the energy spectrum of the fast 
ions. After the interaction, the recoil angle is determined from nuclear 
data, and the recoil energy can then be determined (in two-body reactions) 
by kinematics. Me have added to this calculation the possibility that 

some of the heavy gas nuclei may be concentrated in dust grains. The nuclei 
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considered are those listed in Table II-3 with mean lives longer than 10 
seconds . 

In our calculations we consider spherical grains having an 
exponential size distribution and containing half of the interstellar C, N, 
0, Ne and S and all of the Al, Mg, Si, Ca, and Fe. For such grains, the 
probability that a nuclear interaction takes place in a grain of radius 
between a and a+da is 

VmJLo-- o? / CocC*) 

where a Q is the characteristic grain radius. We assume that the energetic 
particle flux is isotropic, and that the matter density and composition are 
constant throughout the grain. Then if an interaction producing a recoil 
with energy occurs in a grain, we evaluate the grain radius from equation 
II-9, and the position of the interaction in the grain and the recoil angle 



at random. From these parameters, we can then determine the distance r the 
recoiling nucleus must travel to reach the edge of the grain. It is of 
interest to note that the distribution of r in spherical grains can be 
obtained analytically. If the distance from the center of the grain to the 
interaction is x, then the probability normalized to 1 that x lies in dx is 
given by 


O0 < 5 * £ - 32 d-X- (a} 


(ji- id ) 


If we then construct a radius vector through x, the distribution in recoil 
angle 0 about this line is assumed isotropic, or 

"?<***&) du.s = £ itw* 0) 6r-»0 

The cosine law gives r in terms of a, x and cos0 as 

r*- -O.X CC*S & (iT-15.) 

Then the probability of occurrence of r is given by integrating over x„ 


Equation 11-12 implies a minimum value of x exists for a given r, 

X m\n ~ I a ” *"* I (jj - 1 3^ 

Then the desired distribution in r is obtained: 


'?Cr-)d r - rc*) Pi~s») | dr (T- 1'0 

V v, * 

After substituting for x m ^ n » "P(x) and P(cos0), and evaluating d(cos0)/dr 


from equation 11-12, we obtain 
■ a. 


?Crt= J ^ H> 1 (x-,0 

This distribution was given in Ramaty, Kozlovsky and Lingenfelter (1978). 

Once a value for r is determined, we find the stopping range 


X for the excited recoiling nucleus by assuming a Gaussian distribution with 

mean <X(E )> and standard deviation AX as shown in Fig. 11-7, and a lifetime 
r 

from an exponential distribution with mean life given in Table 11-3. If 
the nucleus stops in the grain (r>X) , E * the energy of the recoiling 
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nucleus at the time it deexcites, is evaluated from 


C 


JIB 


t. 


• f> v 
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where p is the density and v the recoil velocity; this case includes the 
possibility that E 7 = 0, i.e. the nucleus stops before it deexcites. 
Equation (19) is also used when r < X, provided that the nucleus , deexcites 




before it reaches the grain edge, i.e,. when 

! Er > t 

0 V dE_ s' L t 
Er 1 j 




Here E 


edge 


, the energy of the nucleus at the edge, is obtained by solving 


- (x(E r }) -r (e -!•?') 

If condition 11-17 is not satisfied, then E * = E e ^S e i.e. the nucleus 

5 x. r * 

leaves the grain in an excited state and deexcites later in the interstellar 

€idg6 

medium where its energy remains essentially E^ ° . The resultant gamma 
ray energy is determined by the velocity corresponding to E^ 7 , and the 
angle between the velocity and the emitted photon, assumed isotropic in 
the nucleus 7 rest frame. 

. The effect of the grains on gamma ray line profiles and the 
dependence of these shapes on the characteristic grain radius are shown in 
Figures II-8 through II-11 0 These figures result from the calculations 
of Ratoaty, Kozlovsky, and LIngenfelter (1978) for an ion spectrum with a 
differential kinetic energy spectral index of 3,’ assuming solar composition 
for the fast ions and the ambient gas. Figure II-8 results from the absence 
of grains; figures II-9, 11-10, and 11-11 show the effect of including 
spherical grains with radii distributed exponentially about characteristic 
fadii of 0.1(1, 0,3(i and 0.5(i, respectively. The case with a = O.lp shows 
very minimal differences from the spectrum with no grains. However, as the* 
characteristic radius increases, there are noticeable enhancements in the 
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very narrow components of the lines at 0.847 and 1.238 MeV from Fe, at 
24 * 2R 

1.36S MeV from Mg, and at 1.779 MeV from Si, as well as in a number of 
lower intensity lines. But the most dramatic and probably easiest to observe 
effect occurs for the 6.129 MeV line of ^0. These figures show that simul- 
taneous observations of a number of these lines with very narrow components 
would then yield information on the size distribution and compositions of 
the grains . 

f) Conclusions 

Using detailed cross sections for charge exchange, ionization, 

and excitation, we have calculated the production rate of Ko; X-rays from 

sub relativistic iron ions. After assuming a set of abundances and energy 

spectra for these ions, we presented X-ray line production rates for both 

a given energy density and a given energy deposition rate of the energetic 

particles. Then, we compared these results with the rate of y-ray production 
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from- the first excited levels of C and Fe. In section d, these calcu- 
lations were applied -to the galactic center, Cen A, and Cas A. For the 
galactic center, the reported observation of the 4.44 MeV line (Haymes et al. 
^1975) are not inconsistent with upper limits on the broad 6.8 keV line. 

However, this consistency with the X-ray data implies that the energy density 

-3 

in fast ions required to produce the observed gamma rays is several keV cm 
.If we accept these results, however, recent observations from HEAO 1 imply 
time variability exists. In the case of Cen A, Hall et al. (1976) reported 
A feature at ~ 4.5 MeV with a width of < 150 keV. The intensity of this 
feature is a factor 0^40 larger than the upper limit on a narrow C line 
set by the X-<ray data. We note, however, that the constraints on gamma-ray 
line intensities imposed by upper limits on broad 6.8 keV emission do not 
apply to sources that are either hotter than a million degrees or opaque 



31 


to the 6,8. keV photons. Thus, the future detection either of nuclear 
•y-ray line emission or of a broad 6„8 keV x-ray line, or both, could give 
very valuable information on the existence of fast ions in the interstellar 
medium and on the physical conditions that lead to the production of these 
particles. In. addition, observations of certain Y“ ra y lines by high resolu- 
tion detectors could provide some constraints on the composition and size 
distributions of interstellar grains. We have shown the dependence of 
these profiles on the grain size distribution in Figs. II-8 through II- 11. 



III. CYCLOTRON FEATURES IN INTENSE MAGNETIC FIELDS 


a) Introduction 

Pulsars with periods in the range from tens of milliseconds to 

seconds are believed to be spinning neutron stars with strong magnetic fields 

rigidly anchored to their surfaces. These values of the rotation period 

are obtained from the conservation of angular momentum during the collapse 

of a main sequence star to a neutron star (with radius ~ 10 fan) . By assuming 

magnetic flux is frozen into the collapsing plasma, the strength of the 

12 

magnetic field is estimated to he ~ 10 G, Pulsations arise from the 
confinement' of the radiating plasma by the magnetic field. Fields of this 
strength alter considerably the interactions of electrons with both radiation 
and ions, primarily by quantizing the allowed electron energies from motion 
transverse to the field. Radiative transitions between energy levels give 
rise to cyclotron line emission and absorption, and the Coulomb scattering 
process becomes highly modified because of the restricted recoil of the 
electron. The line processes have been treated in detail by Daugherty and 
Ventura, (1977 and 1978); the Compton scattering cross section for low energy 
photons (tu « eB/ (m c)) was calculated by Canuto, Lodenquai and Ruderman 
0.971). Resonances in the non-re lativis tic Coulomb electron- ion collision 
frequency were discovered by Ventura (1973). Current X-ray observations, 
however, require more refined calculations of the collisional cross sections. 
In this chapter, we have calculated an effective scattering cross section 
for line photons and the Coulomb cross section for electron-ion collisions, 
treating the electrons relativistically. 

The motivation for the type of calculations carried out in 
this chapter is based on a model for a specific type of X-ray pulsar, a 
highly magnetized star in a binary system with a less evolved companion. 
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X-ray pulsars in binary systems derive their energy from the accretion of 
matter from the nearby companion. The matter is accelerated by the 
gravitational field of the neutron star, and it is believed to be channeled 
by the field onto a magnetic polar cap with zero pitch angle,' There it 
releases a significant amount of its infall energy by heating the atmosphere 
to temperatures at which electrons begin to become relativistic; for example, 
the effective electron temperature of Hercules X—l, the best studied binary 
X-ray pulsar, is greater than 20 keV (Pravdo et al. 1977b). Evidence has been 
presented for an emission line feature at about 60 keV (Triimper et al. 1977 
and 1978, Coe et al. 1977); the electrons involved in producing such a feature 
are mildly relativistic. Thus in this chapter we have considered cyclotron 
line production by electron-ion collisions and the cross section involved in 
effectively scattering the line photons in a relativistically exact treatment. 
In section b, we discuss the quantization of electron orbits in 
strong fields, and show that under the conditions expected at the polar 
caps of X-ray pulsars, electrons will likely spend most of their time in the 
'magnetic ground state, corresponding to zero pitch angle. Then, using a 
fully relativistic quantum electro dynamic (QED) approach, we calculate the 
electron-ion Coulomb cross section and the absorption cross sections and 
emission rates for cyclotron photons. In section c, the Coulomb cross 
sections are used to calculate the energy loss rate of the infalling protons. 
This rate, together with the rate of elastic .nuclear protonrrproton collisions, 
Is then used to calculate the number of line photons an infalling proton 
can be expected to produce, directly or indirectly. We also calculate the 
thermal line production coefficient as a function of the one-dimensional 
temperature for a gas in which the electrons are in the ground state. 

In addition, we find that the infalling protons can create a flux of 





super-thermal electrons streaming down the field lines toward the surface. 
Section d consists of an application of our results to the data from 
Hercules X-l. In particular, we consider the possible observation of the 
cyclotron emission line mentioned above and present a more likely alternative 
interpretation in terms of an absorption feature. Finally, in section e, 
we summarize our results and Indicate tbe need for further observations with 
better energy resolution to distinguish between the models. In addition, 
we propose a more quantitive and comprehensive calculation to determine the 
expected radiation from the whole class of binary X— ray pulsars, 
b) Physical Processes in Intense Magnetic Fields 

The Dirac equation for the spinors -i)r representing electrons in 


an external static and homogeneous magnetic field is, in units where "h = c. = 1, 
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where m is the elect. 


.ectron mass, y^ 1 are the Dirac matrices, and 
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The eigenvalues of total energy W are given by 
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where p Is the electron momentum parallel to the field, c is the speed of 


light, m^ is the mass of the electron, and in. cgs units, 

B* ~~ 4-413 X/0 f3 jau*s. . 


(jx-s) 


The corresponding spinors for a particular positive energy 
gtate, assuming B lies in the z-direction, are 


where h is the normalization length, x ± is the projection of the position 
vector into a plane perpendicular to the field, s is the spin of the electron. 
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and 
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The quantum number j labels the angular momentum parallel to the .field from 
the contributions of spin and orbit around the guiding center, while q gives 
the distance from the guiding center to the z-axis of the chosen coordinate 
system. The functions f^ (x^) are defined most conveniently by the ladder 


operators: 


then 


ir* = x. (■»**; v‘i) 

X* = -fz <.*•**■&)* rg v± '■> 
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define the functions completely. The tt operators correspond to the canonical 
momentum of the. electron and the X operators to the coordinates of the guiding 
center. There are two contributions to the transverse- energy quantum number 
expressed by setting 

^■=.n + 5 + J: , it - Oj ! ; s - - 

Then n refers to the Landau level, the state of the orbital motion of the 
electron about its guiding center, and s is the contribution of spin to 
the magnetic’ moment. Note that only s => -1/2 is permitted for the ground 
state. 


In a magnetic field, radiation can be emitted and absorbed 
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accompanying transitions between transverse energy levels. Under normal 

astrophysical conditions, particle energies are much larger than the energy 

level separation, and the particles consequently have large parallel momenta 

and high quantum numbers. This leads to a classical approximation, and 

while the radiation is still quantized, the large parallel momenta Doppler 

shift the emitted photons into a continuous synchrotron spectrum. However, 
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even for mildly relativistic electrons (y < 1.1) in fields greater than 10 
gauss, the quantum nature of the radiation must he taken into account. 
Discrete cyclotron line emission rates and absorption cross sections have 
been calculated for intense fields by Daugherty and Ventura (1977 and 1978) . 
We have repeated these calculations in the representation above, which is 
somewhat different from theirs. The QED matrix element for a transition is 
given by 


S . = le feL«* f- 1, <x> % c*> 


(UT-iV) ' 


wheye the subscripts f and i refer to final and initial states and (x) 
is the photon vector potential, given by 

f 42fV /l f e ±l 
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where m and k are the photon frequency and wave-vector, s is the polarization 
vector, and V is the normalization volume. The har over i]t indicates the 
Fauli adjoint of the spinor, ijr^ty 0 . 

We have confirmed the results of Daugherty and Ventura (1977 and 

1978) for the absorption cross sections and emission rates using the ladder — 

operator representation in cylindrical coordinates. For example, the cross 
section for an electron in the ground state to absorb a fundamental cyclotron 
photon is given here in the center of parallel momentum frame, where the 
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sum of the parallel momenta of the photon and electron vanishes; 

I ,2m c <i>1 Qj ( [ +C'W’t +M e e^l C£>5 l frfca + {\fj- tfflit 1 ) € a) 
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where 0 is the angle of propagation of the photon and uj its frequency, and 
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6 q and are the polarization amplitudes in the ordinary (wave’s electric 
field parallel to B) and extraordinary ’modes , respectivelyT As we shall 
now discuss., this cross section is of particular Importance for intense 
magnetic fields, since it will determine the opacity in the .line feature. 

In a strong field, mildly relativistic electrons will spend most 
Qf their time in the ground state, unless the densities are very high. To 
see this. Figure III-l shows the lifetime of the first excited state as 
a function of field strength for electrons with no parallel velocity. This 
Was calculated by numerically integrating over all photon transitions to the 
ground state, according to 

% = ( j d 3 K Ko, ct) ) 

•+3 3 

where (k)d k is the rate per electron of photon emission into d k centered 

about k and is calculated in the appendix. As mentioned before, there are 

two possible spin orientations for the excited states; the lifetime of each 

is shown separately in the figure. As can he seen, for low enough densities, 

collisional excitation times will be much longer than radiative lifetimes. 

We will return to this after calculating the collisional cross sections. 


~( 


(Ill'll) 
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In a Coulomb collision, the momentum transfer prefers to occur 
perpendicular to the relative velocity. For energetic protons with small 
pitch angle colliding with an electron in the ground state of an intense 
field, the transfer of energy to the electronic plasma is inhibited (Basko 
and Sunyaev 1975). However, under energetically favorable conditions, the 
proton can transfer transverse momentum to the electron by excitation of the 
quantum levels. In addition, head-on or knock-on collisions can occur, in 
which the electrons are not excited, but are reflected along their field 
line. We have calculated the cross sections for these processes as a 
function of the electron energy in the proton’s rest frame by a single 
.photon exchange approximation. The matrix element is given by 

S+i « -le. JW '* c ho ) <*> i'i «x) (m-ao) 

(c) 

where A q ( x ) is the Coulomb potential due to the ion charge. After a 

Fourier decomposition of this potential as outlined in the appendix, we 


find the cross section <jj for excitation of the jth state from the ground 
state in a Coulomb collision with an ion of charge Ze with zero pitch angle 
is given in the ion’s rest frame by (in cgs units) 

4 J % *>i*Y‘*} £j-(& %fD, Cnt-a) 

where r £ is the classical electron radius, 2.82x10 cm, p and p/ are the 
initial and final momenta to the field respectively, of the electron, E is 
, the electron kinetic energy, and s' is its final spin. There are two 

possible values of the final momentum, given by p' , = + /p a " -2m ^c 2_ B/B T, 

- + — & q 

and these are summed over in equation (3). Finally the £, are defined by 
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The first -three are given here in terms of the exponential integral; 
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The cross section refers to a knock-on collision, where the parallel 

momentum of the electron changes sign (in the ion's rest frame), but no 

transverse excitation occurs. Higher j values refer to collisions in which 

the electron is excited to the j level. These results are shown in Figure 
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III-2 for 2 magnetic field strengths, 10 G and 6x10 G. Because of a 
breakdown in the approximation, the cross sections become infinite at the 
thresholds of excitation. However, the bandwidth of the breakdown is small. 
So we have arbitrarily smoothed these infinities. The approximation breaks 
down because p / , the final momentum, becomes small, and the particles are 
near each other for a long enough time to exchange more than one virtual 
photon. These infinities are the resonances discovered by Ventura (1973). 

As in the non-magnatized case, bound states are even possible. 

From Figure 1II-2, we see that above its threshold, the exci- 
tation to. j « 1 dominates all others, and (p) is roughly inversely 
proportional to B. Thus we can compare the lifetime of the first excited 
state from Figure III-l to the collision time using Figure 111-2 to find 
the density at which electrons may he in an excited state for a significant 


„ time. Explicity the critical density n £x is ...... 

- SXiD 2 * d-m " 3 -for Vsid n cnij-> t (uT'J- 4 ) 

3 

and n increases for increasing B roughly as B . The X-ray emitting regions 
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in binary pulsars are expected to be much less dense than n for B 10 G. 

Under conditions where the electrons occupy the ground state, 


the opacity in the cyclotron line is determined from the cross section given 



in equation III-15. Since this cross section is much larger than that far 
* 

Compton scattering, the possibility arises of a source region optically 
thick to line photons but not to the surrounding continuum. Since Compton 
scattering is the most important factor in radiative transfer outside the 
line process, a line feature will appear in this case. The width of the -line 
can be roughly determined by finding the photon frequency and angle ranges 
where the optical depth exceeds unity. The optical depth is given by 
T b O>> 0) ^ 72e C jr turns’) (d-AS-) 

where d«L is the line of sight increment,. f & (p) is the electron parallel 
momentum distribution. function, and <j is given by equation III-15. The 
asterisk denotes, the center of parallel momentum frame. If we assume that 
the electron spectrum is uniform throughout the source, the 5-function in 
equation III-15 allows us to evaluate the momentum integral to obtain 
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where N is the column density of electrons in the source, p and W are 

the momenta and total energy which satisfy the line resonance condition,. 

given below, and S(u)*, 8*) is the right hand side of equation 111-15 without 

the 6-function. W, and p, are solved for from 
4 - + 

. * \ 
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.which yields: 
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where RHS is the right hand side of equation IXI-27. We note an important 

feature; there exists an angular dependent upper bound w (0) on the line 

max 

energy where the quantity under the square root in equation III-28 becomes 
negative. This limit is given by 

(j4 ^y*j 


(nr -a?) 
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To undertake a study of the radiative transfer of the line 

♦ 

photons, we note that for n n , an electron which is excited by absorbing 
a line photon will deexcite radiatively before anything can happen to change 
its momentum. Thus, the process is effectively a scattering of the line 
photon. The scattering cross section is the .product of the absorption 
cross section and the probability of emission into a given final angle, 
or in cgs units, 

eiCT” (.ul , 6") _ Tf 'g. 4 cJ ' ♦'K u> — \t} { ^ g zvj'i 
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* [(W, Ws+n 1 ^/^ * CVi/ r +vn e c v ) Z .s < rt. 1 3- / £:^ 1 'J J - 
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where everything is as defined for equation 111-15, and the primes indicate 
the emitted photon. The t^’s are shown in Figure III-l . We note the symmetry 
between final and initial photons. The emitted photon has the energy 

■fc tJ ' - [W,- (j +- ^ l/j - J t Cm -3 o 

Solving the expression in the argument of the 6-function in (111-30) for 
0 ) gives the same equation with the primes removed. 

c) Line Production by Collisional Excitation 

Using the results of the previous section for the Coulomb cross 
sections tjj in an intense magnetic field, we have calculated the energy loss 
rate per unit matter traversed for a proton with zero pitch angle by Coulomb 
processes: 

oo * 

jf = ^ ^ <E ‘ flE i CE - 



where m^ is the proton mass, v is the proton velocity, p is the electron 
parallel momentum, f £ is the electron distribution function in p, is 

the relative velocity between the particles in an encounter, and is the 
energy lost in the collision'. We assume" a" relativistic one dimensional 
thermal distribution of electrons, so that f £ (p) is given by 

L <-}>) - K', C { Z4 \ > L" C^wV^/t J (in:***) • 

where K^(x) is the modified Bessel function (Abramowitz and Stegun 1965) and 

and T is the temperature in units of energy. The results are shown in Figure 
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XII-3 for four temperatures and a field strength of 6x10 G. For comparison 

the energy loss rate in a non-magnetized gas of temperature 35 keV is shown 

12 

by the dashed line (Book and Alt 1975) and in a field of 10 G with tempera- 
ture 15 keV by the dash-dotted line. For T = 0, there is a sharp disconti- 
nuity at the threshold for excitation, but this disappears when the electrons 
take on a thermal distribution, and for the higher field strength at 
T = 35 keV, dE/dx is relatively constant in energy, 

Figure III-4 shows the Coulomb mean ranges for the temperatures , 
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discussed above and a field strength of 6x10 G. As Basko and Sunyaev (1975) 

pointed out, for high enough energies elastic proton-proton nuclear scattering 

2 

has a mean free path smaller than the Coulomb range (X , . ~ 50 gm/cm ). 

P~P y C-LclS 

In an elastic scattering, the proton is deflected from zero pitch angle, and 
then the energy can be given up rapidly to ambient electrons, via the excita- 
tion of plasma waves parallel to the field, for example. Thus, the equation 
for the energetic zero pitch angle proton distribution function. in a steady 


state atmosphere can be written 

If 


(pr-Ji) 


Where E is the proton energy, X the amount of matter traversed (in gm/cni for 

example), §■ is the differential (in energy) flux of protons, and cr is the 

PP 
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cross section for a nuclear collision with an ambient proton, taken from 
the compilation of Barashenkov and Maltsev (1961). The solution to this 
equation with the boundary condition that injection occurs at the top of the 
atmosphere, (X=0) with flux cp 6(E-E ) is 


"■i- 38 

To evaluate the line productivity of the protons in a one 

dimensional thermal gas, we have integrated the flux given by equation (III- 35) 
* 

with the cross sections for knock-on and excitation collisions. The production 
rates at a depth X of electrons in quantum level j with momentum p in dp is 
given by ^ 

%i (in-50 

where n g is the electron number density §(E,X) is the proton flux given by 
equation 111-35, v ^(E,p*) is the relative velocity between the particles j 
v(E). is the velocity of the ion, p y * is the incident momentum in the ion rest 
frame, and the momentum in the delta function' is obtained from tbe conservation 
equations ; 

p'* st y (p'rfV') , W'*v WvV/Jf') , 

p± a - = ^'* ( m - m ) 


• ~|W + ) 


where the * T s refer to the ion r s rest frame. The total column density 
production rate is found by integrating along dz, the path of the infalling 
Ions; 


i 


v' 1 J h Vi ^ ’ f 1, ** ?i t f> x 5 


(nr- 3 >) 


The resulting expression for this quantity, after taking the delta functions 
into account, ia 

* (ll-aO 



The excitation of an electron to the jth level was assumed to 

yield j line photons since the electrons tend to deexcite in single step 

transitions (Canuto 1977, Daugherty and Ventura 1977), A knock-on (j=0) was 

assumed to produce as many line photons as' was energically possible, by 

Coulomb collisions with ambient ions. It is important to note that here 

we neglect energy losses of the knock-ons to continuum photons by resonant 

absorption (Pravdo et al„ 1978)-which results in the conversion of these 
/ 

photons to line photons; this process is discussed below. The results for 

the multiplicity of line photons produced under these assumptions are shown 

12 

for the four temperatures with B = 6x10 G in Figure III-5a, and Figure 
III-5b shows the fraction of the proton’s energy put into line photons as 
a function of the injection energy. For Figure III-5b, we assume a line photon 
has an average energy equal to the energy of the 1st excited state with 

p = 0. 

Tn Qrder to investigate the production of cyclotron line photons 
by the high energy tail of a one dimensional thermal electron distribution 
in collisions with ambient protons, we have integrated the cross section 
for the excitation of the first excited state with such, a distribution func- 
tion. The resulting line production rate per unit gas density squared, 
is given by 

■ < (T)* ]4 -f.Cf.) o^lp) trcj.’J. (jr-'ftO 

The results f-or temperatures of 1 to 100 keV are shown in Figure III-6, for 
12 12 

field values of 10 and 6x10 Gauss. As can be seen, the line prqduction 
coefficient rises quite rapidly as kT approaches eB/(m e c) and flattens 
somewhat at higher temperatures. 

As noted above, these results apply to line production by 
electron-ion collisions alone. However, a photon with energy below the 
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cyclotron energy can collide with a fast (P ~0.6 to 0.8) knock-on electron 

and be absorbed as a cyclotron photon because of the Doppler shift. When 

the electron deexcites,- the photon tends to be emitted preferentially into 

the electron’s forward hemisphere, or downward toward the star's surface. 

The importance of this effect may be twofold: reemission at angles around 

90° leads to an increase in line photons, and the absorption results in a 

depletion of the continuum below the line energy, which could make a line 

appear more intense relative to continuum radiation. It is difficult to be 

quantitative about this effect because although the production rate of knock- . 

on can now be calculated, their momentum spectrum depends also on the energy 

loss rate to the radiation, the spectrum of which is unknown. It is, 

however, interesting to note that if the field strength is as high as 

12 

implied by the data of Trtimper et al (1977 and 1978), ~ 6x10 gauss, then 
the knock-on process is the dominant one for producing line photons under 
the previous assumptions (neglect of radiation). So if now we assume that 
the kno f ck-ons lose all their energy by converting ~ 30 keV photons to ~ 60 keV 
photons instead of losing 60 keV per photon produced, then the multiplicity 
is only increased by at most a factor of 2, The knock-on production rates are 

strongly dependent on the field strength because of the threshold effect; at 

\ 

lower field values, relatively more direct excitations than knock-ons occur, 
d) Application to Hercules X-l 

The results of the calculations in the preceding section of the 
fractional energy lost by infalling protons to line photons may be compared 
With recent possible observations of a line feature in the spectrum of 
Hercules X-l at around 60 keV (Trumper et al. 1977 and 1978, Coe et al. 1977). 
Assuming that the line is directly produced by the infalling ions and the 
rest of the ions' energy goes into continuum radiation, the calculations 
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leading to Figure III-5 will apply if radiative', transfer in the source or 

outside does not convert continuum photons into line photons or vice versa. 

There are two observed components of continuum radiation, hard (2 to 20 keV) 

and soft (.1 to 1 keV) X-rays, with roughly comparable luminosities of 1 to 
37 r i 

2xlQ ergs s'” (Giacconi et al. 1973, Schulman et al. 1975, Catura and 

35 -1 

Acton 1975), while the observed line luminosity is ~ 2x10 ergs s . The 

infalling protons are generally assumed to power both continuum components, 

so at most, only 1/2 to 1 percent of the energy comes out in the line. Since 

the effective hard X-ray temperature is 20 keV (Pravdo et al. 1977b), we see 

from Figure III— 5b that this fraction is consistent with infall proton energies 

of 140 MeV to 180 MeV. However, we still must investigate the effect of 

radiative transfer, in the source and outside. 

To ascertain whether radiative transfer is important, we need to 

estimate the electron density in the source. From Figure II 1-4 , we see that 

-2 

the range of accreting matter is at least of order 10 gm cm , corresponding 
24 - -2 

to «■- 6x10 a cm . The reported width of the line feature is 12 keV (FWHM) , 
and this places an upper limit on the linear dimension of the source parallel 


to the field, since the line energy varies as field strength. The upper 

_3 

.limit is obtained from assuming B “ R as for a dipole field: 

L-s i * 1r - 3* tr * 7x/a v ctn (in -it) 

where R is the neutron star radius. Thus the lower limit on the average 


(M ' Vi) 


electron, density in the source is given by " 

From equations 111-25 and 111-15, then, we estimate the opacity to be 

(airm e t v T [lokm) (-if kel) V® J 
This implies that if the line of sight covers 50cm or more of the source. 


radiative transfer will be important 
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The preceding considerations apply if the infalling protons 
directly excite the line. If, however, the protons simply thermalize without 
yielding appreciable line emission, the problem becomes that of calculating 
a thermal production coefficient. Applying the results of such a calculation 
shown in Fig. 1II-6 to the observation of Trumper et al. ,(1977 and 1978), 
we obtain the product of the. density squared and the volume V of the source 
region 

=• SLX/D S1 ’ c™* 3 . 

12 

for a temperature of 35 keV and field strength of 6x10 G. As before, we 

.can set an upper limit on the volume, and thus a lower limit on the density, 

-by -assuming that the_ transverse area occupied by the accreting matter is at 
2 

mpst tt • (0.1R) , or 

' A <- 3x/d"W Clfc.j' ' 

.This limit, together with that expressed by III-38, implies 

■n * 3x/o“W (jfz,) 3 ' Car-vO 

implying that transfer is important in this case, also. Therefore, an 
optically thin treatment of this radiation, as in Daugherty and Ventura 
(1977) is not likely to relate directly to observations. 

have considered the problem of line transfer in a rudimentary 
way, by calculating the opacity at a given angle and frequency for a thermal 
gas of electrons. The motivation lor this is the following: where the 

opacity is large, photons are trapped in the source and must diffuse out, 
qnd at the surface, the intensity radiating outward approaches approximately 
its blackbody value. Therefore, where t(id, 9) 1, we assume the line rises 

out of the continuum. The upper limit to the width of the line reported 
by Trumper et al. (1977 and 1978) is 12 keV; as-mentioned above, the field 
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Implied is ~ 6x10 G, and in addition, they show evidence that this feature 
is produced in a gas with temperature ~ 35 keV. So for these conditions 
and four values of electron column density, we have estimated the width of 
a cyclotron line feature for several viewing angles and column densities by 
finding the frequencies where t - 1. Since Doppler shifting by the motion 
of the' electrons is responsible for most of the smearing of the opacity, we 
expect the line to he narrowest at 90° to the field lines (and electron 
motion) and to broaden as the angle departs from this. Another - factor 
which is important is the maximum frequency allowed, given by equation III-29, 
especially near 90°, where it is responsible for the high energy drop in line 
opacity. 

The results for the line width are shown in Table III— 1, for 

22 23 

viewing angles from 10° to 90°, and electron column densities 10 ,10 , 

24 25 —2 

10 , and 10 cm . As can be seen, for an electron temperature of 35 keV, 

this calculation cannot reproduce the width of the observed feature. The 

22 -2 

column density can’t be much lower than 10 cm because of the estimates 
Of lower limits on the density expressed in equations 111-39 and III-43. Thus 
to be consistent with the observation, we are forced to assume a lower 
electron temperature. As mentioned before, lower energy measurements have 
placed the temperature greater than 20 keV; so we have shown the line width 
results for this value also. Here we find that the reported width can be 
reproduced if the viewing angle is restricted to be very close to 90° and 
if the column density is the very minimum allowed. We note that the same 
conclusion about the yiewing angle has been arrived at classically (Meszaros 
1978). We also point out that the line of sight through the 20 keV source 
rjust be less than 100 cm, at most. 

There is an additional problem with the interpretation of the 
observation as a cyclotron line feature because of its pulse phase dependence. 



The 1.24 second pulse period shows about a factor of 4 increase in 2 to -20 
keV X-rays over about 30% of its phase (Pravdo 1976), Because of the sudden 
rise and fall of the intensity pulse, the most likely explanation is occulta- 
tion by a corotating cloud of gas which is presumably responsible for the 
soft X-ray flux mentioned above. The luminosity and temperature of this soft 

g 

source imply that it has linear extent, of around 10 cm, which places it well 
above the surface of the star. A likely candidate is a corotating Alfven 
shell, plasma concentrated in the region where the magnetic field begins to 
dominate the inflow (McCray and Lamb 1976, Basko and Sunyaev 1976). The 
observations of Trumper et al. (1977 and 1978) together with the analysis of 
Kendziorra et al (1978) imply that the line is pulsed in phase with the 
'2' to 20 keV k-rays. Since we 'have shown that observation of the line at 
angles smaller than ~ 85° will contribute significant broadening (see Table 
1II-1), this implies that the pulse is a very narrow fan beam. If this were 
the case, then the Alfven shell would have to suddenly thin at an angle 
close to 90® to the polar field lines, the opposite of what has been proposed 
(McCray and Lamb 1976, Basko and Sunyaev 1976). Thus we conclude that this 
interpretation is unlikely. 

A more likely interpretation of the "feature 1 ' at 60 keV which is 
still consistent with the data is that it is simply residual continuum 
radiation with an accompanying absorption dip at lower energies. The absorp- 
tion would be due to the magnetic field, but the assignment of a field 
strength may not be straightforward because of Doppler shifts. No such 
shift would result if the absorbing electrons were cold, but because of the 
high field strength required, they must be located in or near the source; 
consequently, they are bathed in X-rays and should approach an equilibrium 
temperature. In' this case, they become part of the blackbody, and no 
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absorption will appear. The knock-ons discussed in the preceding section 
are a possible solution to this problem. These electrons stream down the 
field lines with P ~ 0.6 to 0.8, and can consequently absorb upward moving 
photons below the cyclotron frequency to as low as 20 keV. Since they are 
mildly relativistic, they preferentially reemit a line photon in their 
forward directibn, i.e. downward. If the column density of knock-ons is 
such that the absorption depth doesn’t exceed unity, then re-emission into 
the backward hemisphere will be not important. Consequently, a depletion of 
photons below the cyclotron frequency results, appearing as an absorption 
dip, but the field strength is not readily obtainable until the electron 
spectrum is known. 

This picture predicts a certain behavior of the absorption feature 
with pulse phase. In particular, as the line of sight comes closer to the 
polar field lines, the energy at which absorption begins should become lower, 
and the variation in energy should become smoother. In a recent paper 
(Pravdo et al. 1978), we assumed an arbitrary high energy electron distribu- 
tion, and directly compared data from the cosmic X-ray experiment on OSO-8 
**ith these predictions. As discussed there, the fit was acceptable over the 
region of the spectral pulse, and the best fit angles generally performed 
20° excursion in a region around 40°. As opposed to the line emission 
Interpretation, however, here the spectral pulse (which may only be loosely 
correlated with the Intensity pulse probably formed by the Alfven shell) 

defines a pencil beam (or cone). The best fit field strength was coinci- 

12 

dentally consistent with the emission measurements, ~ 6x10 G. However, 
the cyclotron photon energies here are too' broad to form a well-defined line 

e) Summary and Discussion 

12 

The strong magnetic fields (> 10 G) believed to exist at the 


feature. 
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• surfaces of neutron stars can strongly modify the interactions of electrons 
with ions and radiation. We have considered the effects of such modifications 
on the X-ray spectra of the pulsars in binary systems which derive their 
luminosity from an accretion from their companions. The effect of the 
intense field is to quantize the allowed electron momenta transverse to the 
field lines. Radiative transitions between levels give rise to cyclotron 
emission and absorption, and since transverse momentum transfer from an ion 
to an electron must occur in discrete quanta, the Coulomb cross section becomes 
a knock-on cross section and a series of excitation cross sections. We have 
shown that electrons in the polar caps should almost always be found in their 
ground state and have then calculated the cross sections for. these processes 
treating the electrons relativistically . 

Since Hercules X-l is the best studied of these binary X-ray 
pulsars and shows a possible line feature, we have used the above cross 
sections in an attempt to understand the physical conditions in the polar 
cap of this source. Data from 2 to ~ 20 keV imply that the effective electron 
temperature in the source is greater than 20 keV, To produce the observed 
line feature with such a narrow width, we find that the temperature cannot 
be higher than 20 keV, the angle between our line of sight and the polar cap 
field lines must lie within 10° of 90° during the entire high intensity part 
of the pulse, and the linear extent of the source along the line of sight must 
be less than 100cm, Since these requirements are quite stringent, we have 
attempted an alternate explanation of the data. Due to the limited resolution 
of scintillation detectors, the observation could not distinguish between 
an emission line and an absorption dip at lower energies, although the authors 
preferred to stress the former interpretation. The absorption dip could be 
more likely, however, hecause of the ability of the ions to produce a high 
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parallel momentum population of knock-ons, which would absorb upward moving 
continuum photons and reemit them preferentially downward. This interpre- 
tation has the added feature that it provides an acceptable fit -to- the observed 
steepening at ~ 20 keV. 

An observation of this source with high resolution detectors 
would discriminate between these possibilities. We note that some clarification 
should be forthcoming from a recent point at Her X-l by the HEAO-1 A2 experi- 
ment. In addition, there are several other objects in this category and 

s 

presumably we are seeing them at different ranges of viewing angles. High 
resolution spectroscopy of many of them can probably be used to gain a 
comprehensive understanding of the physical processes occurring in the 
strong magnetic fields of these objects. 



IV. THE ■ ANNIHILATION OF GALACTIC POSITRONS 
a) Introduction. 

Since positrons from outside the heliosphere have been observed 
at energies from about 200 MeV to 10 GeV, and a gamma ray line at 0.51 MeV 

has been observed from the galactic center direction, it is reasonable to 

' • * 

expect that there may be a significant number of positrons annihilating 
in the galaxy. The primary sources" of these positrons are the decay of 
nt mesons, 3*" decay of radioactive nuclei and perhaps the strong electric 
and magnetic fields of pulsars. The tt + mesons are created in energetic 
particle collisions; the threshold in a proton-proton collision is ~ 290 MeV 
if one proton is at rest. Most of, the time, this meson decays into a p + . 
and its neutrino; the p"*" then decays into a positron. Positrons produced in 
this way by galactic cosmic rays probably account for a significant fraction, 
if not all, of those observed at earth. The radioactive nuclei which decay 
by positron emission can be formed in nuclear reactions involving protons 
and alpha particles on C, N, and 0, as tabulated by Ramaty, Kozlovsky, and 
LIngenfelter (1975). These reactions are expected to occur in regions with 
large fluxes of low energy cosmic rays, which may be concentrated around 
supernova remnants. Also, heavier decay nuclei with long lifetimes can 
be formed in the explosive nucleosynthesis associated with supemovae and 
noyae,' In either case, supernova remnants are a likely site of positron 
concentrations . The two production modes can possibly be distinguished, 

since low energy cosmic rays produce a strong 4.443 MeV line, while explosive 

' 26 22 26 
nucleosynthesis produces A1 in supemovae and Na in novae. The A1 

26 22 22 

decays to Mg and the Na to ^Ne excited states, and they deexcite by 
emitting gamma rays at 1.809 and 1.275 MeV, respectively. If any of these 

;w. i i 

lines were observed accompanying the 0.511 MeV line, the production mode of 
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the positrons could probably be determined. Finally, pulsars have been 

suggested as a positron source (Sturrock 1971), In their intense magnetic 

fields, gamma rays of energy greater than 1.022 MeV will quickly produce 
+ — * 

e -e pairs, and if a strong electric field is present, the positrons can be 
accelerated away from the pulsar before annihilating. All these sources 
• produce relativistic positrons. 

Leventhal, McCallum and Stang (1978a, b) have recently reported 
observation of positron annihilation radiation from the Galactic Center using 
a balloon borne germanium detector. The observed line is at 510.7 +0.5 keV, 
and its full width at half maximum (FWHM) is less than 3.2 keV, The total 
flux in the line was (1.22 +0.22)xl0 ^ photons cnf^s There is also some 
evidence for the three-photon continuum from triplet positronium annihilation. 
The 0.511 MeV line was previously seen from the solar flares of 1972, August 
4 and 7 (Chupp, Forrest and Suri 1975), but because the lower energy resolu- 
tion of the Nal detector used, only an upper limit of about 40 keV could be 
set on the line width from this observation. 

Depending on the temperature and density of. the ambient medium, 
positrons and electrons can either annihilate directly or form positronium. 

The importance of positronium formation in the interstellar medium was pointed 
QUt by Steigman (1968), by Stecker (1969), and by Leventhal (1973), and 
positron annihilation in solar flares, both direct and via positronium, has 
been treated by Crannell et al. (1976). Positronium can fort* either in the 
singlet state which annihilates into two 0.511 MeV photons, or in the triplet 
state which decays by 3 photon annihilation. 

Positronium is formed by both radiative recombination with free 
electrons and charge exchange with neutral hydrogen (atomic or molecular) , • 
Charge exchange with heavier ions is much less important than these processes 
(Crannell et al. 1976), • Once formed in a particular spin state, positronium 
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in the interstellar medium annihilates from the same state, because -the 
lifetimes of both singlet and triplet positronium (10 ^ sec, and 10"^ sec, 
respectively) are much shorter than typical collision times with interstellar 
gas , 

Except for ultrarelativxstic positrons, the probability for 
annihilation in flight is negligible until the positrons slow down to 
energies of the order of several hundred eV. (Ultrarelativistic positrons 
have a small probability ~ 10%, for annihilating in flight, but these 
annihilations produce only a broad continuum which in general cannot be' 
observed above other continua.) The purpose of the present paper is to 
investigate the annihilation process at energies below a few hundred eV. 

At these energies positrons either form positronium in flight or thermalize. 
The importance of free electrons in thermalizing positrons was first pointed 
out by Crannell et al. (1976). 

In Section b we discuss the energy loss processes of positrons 
in a partially ionized medium, and the charge exchange cross section with 
neutral hydrogen; in Section c we calculate, using. a Monte Carlo simulation, 
the probability of positronium formation in flight in a partially ionized 

I 

gas as a function of temperature, density and degree of ionization; in 
Section d we present the rates of direct annihilation and positronium forma- 
tion of thermal positrons as a function of the temperature; in Section e 
we calculate line shapes resulting from positronium formation in flight and 
after thermal ization; in Section f we discuss the implications of the 
observed line width on the site of positron annihilation in the Galaxy and 
the consistency of the flux in the line with various positron production 
modes; and we summarize our results- in Section g. 
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b) Energy Losses and Charge Exchange of Positrons 

When a positron is injected into a partially ionized medium, it 
loses energy to. the. ionized, .component by exc-i-ting plasma - waves and 'to the' 
neutral component by the excitation and ionization of hydrogen atoms and 
molecules. The energy loss rate in the plasma, measured in eV/cm, is given 
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by (Book and Ali 1975) 

= i.SxiD ' 13 nc. - M'CtOl 

Where E is the positron energy in eV, is the electron density in cm , 
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T is the temperature of the electrons in the plasma, 
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and M* is the derivative of M with respect to its argument. 

Unlike the energy loss process in a plasma, which is essentially 
continuous, in a neutral medium the positron loses a significant fraction of 
its energy whenever it excites or ionizes an atom or molecule. The cross 
sections for ionizing ‘atomic hydrogen- and exciting its 2p level from the 
ground state by electrons are shown by solid lines in Figure IV-1 (Kieffer 
1969). The cross sections for exciting other levels are small in comparison 
to that for the 2p level (Kieffer 1969, Moiseiwitsch and Smith 1978), and 
are not expected to contribute significantly to the energy loss. We assume 
that the cross sections of these processes for positrons are the same as 
for electrons. Also shown by a solid line in Figure IV-1 is the cross 
section for positronium formation by charge exchange with atomic hydrogen. 
This cross section was calculated by Drachman, Omidvar and McGuire (1978) 
for positronium formation in the ground state. The inclusion of the excited 
states may increase this cross section. 
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The dashed line in Figure IV-1 are positron cross sections in 
molecular hydrogen. The total cross section is from laboratory measurements 
made by Kauppila et al. (1977). (Similar measurements by Coleman, Griffith, 
and Heyland (1974) give somewhat lower values). The ionization cross-section 
is from the compilation of Keiffer and Dunn (1966) for electrons in molecular 
H, for the charge exchange cross section above 50 eV, we use the calculation 
of Sural and Mukherjee (1970). This cross section would also be significantly 
larger if the excited states make an important contribution. Below the 
ionization and excitation thresholds we have estimated the charge exchange 
cross section by subtracting from the total the positron elastic scattering 
cross section of Hara (1974). We then connected the low energy values with 
those above 50 eV as shorn in Figure IV-1. We assume that above the excita- 
tion threshold, the difference between the total cross section and the sum 
Of those of ionization and charge exchange is the excitation and dissociation 
cross section of moleculdr hydrogen by positrons. We wish to note that except 
for the total, the cross sections in molecular H are quite uncertain and we 
call for further studies to determine their values. 

In atomic hydrogen, the energy lost in an excitation collision 
is 10.2 eV. For ionization, the mean energy of the ejected electron is 
1/4, of its binding energy, and the distribution of the ejected electron 
•Velocities can be approximated by a Gaussian with standard deviation 
0,382 ac ( a *» 1/137) (Omidvar 1965). In molecular hydrogen we assume an 
energy loss of 12 eV in an excitation and dissociation collision (Herzberg 
1950) and we take the same distribution of the ejected electrons as in 
atomic hydrogen except for the different binding energy, 
c) Positronium Formation in Flight 

We consider a positron of energy E q injected into a partially 
ionized hydrogen gas of temperature T and ionization fraction X = n^/ (n e +n^) , 
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where and are the electron and neutral hydrogen densities. If the 
temperature of the ambient medium is of order 10^K or larger, the positrons 
will thermalize before forming positronium. If, however, the temperature is 
lower, then some fraction of the positrons form positronium in flight by 
charge exchange with neutral hydrogen. To evaluate this fraction, we have 
performed a Monte Carlo simulation in atomic hydrogen with ionization 
fraction X and T < 10 K. Slowing down of positrons in molecular hydrogen is 
treated' separately below. 

We then allow the positron to lose energy continuously to the 
ionized component until it makes an inelastic collision with a hydrogen atom. 
The energy, E^, at which this collision takes place is determined from 

71 

where R is a uniformly distributed random number, dE/d-L is giyen by equations 
(iy-1) through (IV-4), and a is the sum of the ionization, excitation and 
charge exchange cross sections in atomic hydrogen shown in Figure IV-1. The 
probability that any one of these processes occurs at is proportional to 
its cross section at this energy. If positronium is formed, we remove the 
positron from the calculation because positronium in both the singlet and 
triplet states is expected to annihilate in a time much shorter than its 
collision time in the interstellar medium. In an excitation or ionization 
collision on the other hand, the positron loses energy as discussed in 
Section b. After an excitation or ionization we repeat the process specified 
by equation Ciy*'5) with a new value for E q given by the difference between 
and the energy lost in the collision, and we continue this procedure 
until either positronium is formed or the energy falls below the positronium 
formation threshold. 
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The fraction of positrons which form positronium in flight before 

thermalizing, f , is shown in Figure IV-2 as a function of X for two cases; 
ps 

n e = 0.1cm ■*, T = 8000K and “ SxlO^cm \ T = 1.16xl0^K. The parameters 

of the first case could be representative of the warm component of the 

interstellar medium (McKee and Ostriker 1977), while those of the second 

case are essentially the same as the values considered by Crannell et al . 

0,976) for solar flares. The fraction of positrons forming positroniuni 

in flight increases with decreasing ionization fraction; for the interstellar 

case and X < 0.05 more than half the positrons form positronium before they 

•thermalize. If n £ = 0, about 95% of the positrons form positronium in 

flight. We also see that f increases with increasing electron density 

P s 

because of the reduction of the energy loss rate to the plasma due to the 
reduction of the Coulomb logarithm, whose argument is given in equation (IV-3) . 
We have calculated f for temperatures other than shown in Figure IV-2, 
and 'we £ind that it varies by less than 20% for T < 6xl0^K. For higher 
temperatures, f ^ is expected. to be small in any case because of the absence 
of neutral hydrogen. 

We have compared the results of Figure 2 with the calculations 
of Crannell et al. (1976) and we find that for X = 0.5, our f is smaller 

ps 

by about a factor of 2 than theirs (see Figure 4 in Crannell et al, 1976). 

This results from an erroneous factor of 2 in the Fokker-Planck equation 
that was used by them to treat the slowing down of positrons in a partially 
ionized plasma (G. Joyce, private communication 1978). For X = 0,09, the 
discrepancy is much smaller because in this case, energy loss to the plasma 
is less important relative to charge exchange. We wish to note that for 
low values of X, the effect of the neutrals on the energy loss becomes 
important; for example, for X n 0.09, neglecting the neutrals would Increase 
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f by about 20%, 
ps 

We do not expect molecular hydrogen in a phase of the interstellar 

medium with high temperatures and. high ionization fraction. -For molecular 

hydrogen, therefore, we only evaluate f for - 0. Using the same Monte 

Carlo simulation as for atomic hydrogen and the cross sections of Figure 

iy-1, we obtain" f “ 0.93," This yalue is slightly lower than f in atomic 
' ps- - w ps 

hydrogen at n^ = 0 because of the higher threshold for positronium formation 
in the molecular case. We note that the possible increase of the charge 
exchange cross sections due to the inclusion of the excited states would 
increase f for both atomic and molecular hydrogen, 
d) The Fate of Thermal Positrons 

Having evaluated the fraction of positrons that form positronium 
in flight, we now consider the fate of those positrons which thermalize 
before annihilation. Thermal positrons will eventually form positronium 
or annihilate directly, and both these processes can occur with either hound 
or free electrons, Positronium formation with free electrons is radiative 
recombination. Radiative recombination and direct annihilation rates with 
free electrons were calculated by Crannell et al. (1976); we have extended 
these calculations to lower temperatures, and the results are shown in 
figure IV-3. The direct annihilation rates with atomic .electrons were 
Obtained by Bhatia, Drachman and Temkin (1977) from a detailed quantum- 
mechanical scattering treatment. Their thermally averaged annihilation 
rates are also shown in Figure 3. The rate of positronium formation due to 
radiative recombination includes the contributions of the excited states 
(Seaton 1959), 

The rate of charge exchange of thermal positrons with residual 
neutral hydrogen can be calculated by assuming that the positrons are in 
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thermal equilibrium with the ambient medium. This is justified because even 
though annihilation and positronium formation tend to modify the positron 
distribution, the rate of reestablishing thermal equilibrium is much faster 
than the sum of the rates of all annihilation processes. Therefore, the rate 
of charge exchange of thermal positrons can be written as 

DO 

TW (E) CQ ce') tTCE'i dE. (ut-O 

where n^ is the Maxwell-Boltzmanh distribution, and a ce is the charge 

exchange cross section with atomic hydrogen shown in Figure IV-1, We have 

evaluated equation IV-6, and the results are shown in Figure TV-3. As can 

be seen, this rate is a very strong function of temperature, but even at 

its maximum, it is lower than the thermal ization rate (Spitzer 1962) . 

'From Figure IV— 3 we see that for temperatures greater than 10^K, 

the dominant annihilation process is direct annihilation with free electrons . 

At lower temperatures, radiative recombination becomes more important than 

direct annihilation, hut depending on the amount of residual neutral 

hydrogen present, charge exchange may in fact be the most important process. 

5 

JPor example, at T = 10 K, the rate of charge exchange is greater than that 

—6 ■ 

of radiative recombination if n^/n^ > 5x10 . Since for temperatures less 

5 

than 10 K, we expect larger neutral hydrogen concentrations, charge exchange 

4 

should be the dominant process down to temperatures around 10 K. At 8000 K, 
R ce /n^ - R /n g , so that if, for example, the medium is half ionized (e.g. 
the warm ionized component of the interstellar medium, McKee and Ostlriker 
1977), the rates of charge exchange and radiative recombination are about 
equal. , At lower temperatures, radiative recombination should remain the 
dominant process as long as there is an appreciable concentration of free 
electrons. However if n /n becomes small, a large fraction of the positrons 
form positronium in flight, as can be seen from Figure IV-2. 
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From Figure III-3 we also see that there is a finite probability 
for direct annihilation even when positronium formation is the dominant 
process. For example, at T = 8000K and n £ = n^, about 7% of the positrons 
will annihilate directly (mostly with free electrons). 

e) The Energy Spectrum of Positron Annihilation Radiation 

In this section, we evaluate the energy spectrum of positron 
annihilation radiation resulting from annihilation in a cold neutral and a 
partially ionized medium. We first consider the cold neutral medium where 
the radiation is produced by positrons forming positronium in flight and by 
positron's annihilating directly with bound electrons. Using the Monte Carlo 
singulation described in Section c, we calculate the- photon spectrum resulting 
fronj positronium formation in molecular and atomic hydrogen with = 0. 

■The probabilities for forming positronium in the singlet or triplet state 
are 0.25 and 0,75, respectively. Annihilation from the singlet state produces 
two photons of 0.511 MeV in the rest frame of the positronium. In the frame 


of the galaxy, the energy of a photon from such an annihilation is 

I* 




Cffl-7) 


where £ = 0.511 MeV, P is the velocity of the positronium, and 0* is the angle 
between the positronium velocity and the photon direction of the positronium 
rest frame. We assume that the photon direction is isotropic in this frame. 
Annihilation from the triplet state produces three photons whose energies 
can also be obtained from equation (IV-7) proyided that is replaced by 
•£* which is chosen according to the distribution given by Ore and Powell 
(1949). The results for positronium formation are shown in panels a and b of 
Figure 1V-4, and as can be seen there is little difference between the line 
shapes fot the atomic and molecular cases. The full width at half maximum 
for both these cases is about 6.5 k^V. 
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For direct annihilation with atoms or molecules, the line shape 
is determined by the momentum distribution of the bound electron modified 
to account for' the fact that the positron cannot approach the nucleus too 
closely. For the positron-hydrogen atom system, this momentum distribution 
was accurately calculated by Rumberston and Wallace (1972), and from it the 
gamma-ray line shape is easily obtained. The FWHM is only 2.6 keV, consider- 
ably less than that from positronium formation in flight. For the positron- 
hydrogen molecule system, there exists an experimental determination of the 
momentum distribution (Briscoe, Choi and Stewart 1968), and a FWHM of 3.2 keV 
can be extracted from their work; however, since their measurement was done 
in liquid H^, it is expected to be slightly too wide. The reason for this 
is discussed by Humherston (1974) in the case of liquid helium. The scaling 
of his results to molecular hydrogen results in an expected FWHM of only 
3 keV; consequently we have reduced the energy scale of the work of Briscoe 
et al. (1968)- by a factor of 3/3.2. 

Panels c and d of Figure 1V-4 show the spectrum of positron 
annihilation radiation in a cold neutral medium. The direct' annihilation 
spectra were- integrated. over 1 keV intervals and added to the positronium 
formation histograms shown in panels a and b. Although only 6,5% and 5.5% 

Of the positrons in molecular and atomic hydrogen, respectively, fall below 
the Ps formation threshold (from the Monte Carlo simulation described above) , 
their contribution is enhanced because direct annihilation always yields 
two gamma— rays , while Ps decay yields three photons 75% of the time. -As can 

I 

be seen in Figure IV-4, the FWHM of the 0.511 MeV line in cold neutral 
hydrogen is about 5 keV. If, however, the inclusion of excited states 
increases significantly the charge. exchange cross section, then fewer 
positrons fall below the Ps formation threshold. In this case, the line 
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width is closer to the value (6.5 keV) determined for positronium formation in 
flight. The effect of the three photon annihilation is manifest in a slight 
asymmetry in the line shape ^nd a low energy tail , 

In the case of annihilation in a partially ionized medium 
(X > 0.05), most of the positrons thermalize rather than charge exchange in 
flight. As discussed in the previous section, thermal positrons will directly 
annihilate and radiatively recombine with free electrons, and if the gas is 
warm enough, they form positronium by charge exchange with residual neutral 
hydrogen. The shape of the 0.511 MeV line from both singlet positronium 
annihilation following radiative recombination and direct annihilation 
with free electrons is expected to be a Gaussian of FWHM =- 0.011 (keV) 

(T(K)) ' (Crannell et al. 1976). For example, if T = 10 K, the width 
is 1.1 keV. 


The photon spectra, q (e ) , resulting from charge exchange of 
thermal positrons was calculated as follows: For annihilation from the 

singlet state 


/f S ( g^CE) TT dt. 

2 

where =* + 6. 8 eV, while for the triplet, 

% (£,), [“je . P**,. 

*P ‘ 1 


x 


(jz .-?) 


(&-■>) 


Where 


E c rE2^ or £ Y > £ ° 

2 l 6T8eV for £ e 0 * 

For equation (8), P t (x) is the distribution of Ore and Poweli (1949) with 
X = e */e e> Xj, = £y/ (s 0 Cl+P)) , and = min (1, £^/ (ye o (1-8) ) , where y 
and 8 are the Lorentz factor and velocity (divided by c) of the positronium 
atom formed, corresponding to (E -6,8 eV) , 

By combining the results of equations (IV-8) and (IV-9), we have 
plotted in Figure IV-5 photon spectra from thermal positron annihilation via 
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positronium formed by charge exchange for T = 8000K. and 5000K. The FWHM’s 
-for these two cases are about 1.5 and 3,2 keV, respectively, 
f) Discussion 

As we have pointed out in the Introduction, the FWHM of the 0.511 
MeV line from the Galactic Center was observed to be less than 3.2 keV 
(Leventhal et al> 1978a, b) . When compared with the line shapes shown in 
Figure IV-4, this width is narrower than that resulting from annihilation 
in a cold medium. However, in a partially ionized medium (X p 0.05), for 
temperatures less than ~ 10^K the calculated line width is less than the 
upper limit set by the observations. This result implies that the medium 
in wMch the positrons annihilate should be substantially ionized. 

If we consider the three component model of the interstellar 
medium (McKee and Ostriker 1977), then the most acceptable sites of annihil- 
ation are the warm neutral and the warm ionized media, whose. temperatures 
are about 8000K and ionization fractions are 0.15 and 0.68, respectively. 

For this temperature, the width of the line is only ~ 1.5_keV, as can be 
seen from Figure IV-2, 80% and 98% of the positrons thermalize for these 
choices of ionization fractions. The bulk of the positrons probably do not 

annihilate in cold clouds even though most of the interstellar matter is 

-5 

believed to he in them. Since n e /n^ is very low in clouds 10 , Guelin 

et al. 1977), the line width would probably be larger than observed. Neither 

should they annihilate in the hot component whose temperature is about 

5x10 K and density 3x10 cm ; however, we do not expect much annihilation 

in this medium in any case, because the annihilation time is on the order 
8 

of 3x10 years, and in this time the positrons are expected to encounter 
the denser clouds in which their annihilation time is much shorter. 
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Other possible annihilation sites are young supernova remnants, 

* 4 

For example in the filaments of the Crab Nebula, n £ =- and T £ ^ 10 K 

(Davidson and Tucker 1971). Supemovae are also likely sources of positrons 

through all the. processes' mentioned in the Introduction. 

The 0,511 MeV line was observed from the direction of the Galactic 

Center, but because of the wide opening angle of the detector (~ 15°) the 

.field of view included a large portion of the Galactic Center region extending 

to the expanding arm of 135 km/sec (Mezger 1974). If the positron source, 

however, is located at the nucleus of the Galaxy, the positrons would 

therraalize and annihilate in the extended HII region in the nuclear disk, 

✓ * 

even though most of the diffuse- matter is in the 270 pc molecular ring 

(Scoville 1972) encircling this- region (Mezer 1974). The ionized region has 

-3 4 

an electron density of 15 to 30 cm and temperature less than- 10 K. From 

equations (IV-1) through (IV~4), we expect positrons with initial energies 

of seyeral hundred keV to traverse a total distance of about 500 pc in about 

2000 years before they thermalize. Because of scattering, however,, the 

rectilinear distance could be much smaller, and therefore ' the positrons 

would annihilate before reaching the molecular ring. The annihilation time 

3 

in the ionized .medium is about 10 years. For the temperature of this medium, 
the width of the line is less than about 1.5 keV, consistent with the observed 
Upper limit of 3.2 keV. 
g) Summary 

We have investigated the annihilation of galactic positrons, and 
the formation of the 0.511 MeV line and the accompanying 3-photon positronium 
continuum. In a ,cold and neutral medium, about 95% of the positrons form 

positronium atoms in flight and these atoms annihilate before undergoing 

- . 15 _3 

further collisions if the density of the medium is less than ~ 10 cm 



The shape of the 0,511 MeV line from positron annihilation in a cold neutral 

medium has a FWHM of about .5 keV (Figure IV-4 panels c and d) , and this 

width appears to be inconsistent with the observed width of the 0,511 MeV 

line from the Galactic Center (< 3,2 keV), Leventhal et al. 1978a, b)'. 

In a warm, partially ionized gas of temperature less than several 

times 10 K, only a .fraction of the positrons form positronium in flight and 

this fraction decreases with increasing n^/n^ (Figure IV-2). At higher 

temperatures, essentially all the positrons thermalize before forming 

positronium or annihilating directly. 

The rates for the various processes leading to thermal positron 

annihilation are shown in Figure IV-3. For temperatures less than ~ 10^K, 

positronium formation is the dominant annihilation channel, while for 

.higher temperatures free annihilation dominates. The width of the 0.511 MeV 

line from thermal positron annihilation is about 3 keV for T =“ 10** k and 

A 

yaries as approximately T 2 . The observed .width from Galactic Center suggests 
a temperature less than this value, and a degree of ionization larger than 

i 

about 5%, for which half the positrons thermalize. For these conditions more 

i 

than 90% of the positrons annihilate after forming positronium, and this is 
consistent with the possible evidence for 3-photon continuum found by 
leventhal et al. (1978a, b). 

Among the various galactic sites, the warm component of the 
interstellar medium, filaments or knots in supernova remnants, and the 
extended HII region in, the nuclear disk could be sites for forming the 
0,511 MeV line from positron annihilation. 



V. SUMMARY 


This work has been concerned with emission lines and other spectral 
features, in the hard X-ray and soft gamma ray spectra of various astrophysical 
sites. In this chapter, we summarize the results and attempt to point out 
the value of future calculations and experiments which could yield information 
on such sources. Many of the observations could be carried out from the 
proposed Gamma Ray Observatory (GRO), a satellite which would carry a comple- 
ment of gamma ray detectors including one of high spectral resolution. 

In Chapter II, we have calculated in detail the expected emission of 
a broad (FWHM «2.4 keV) line at 6.8 keV from the 2p to Is transition of low 
energy iron ions. It was shown that this line is sensitive to ions at 
~ 10 MeV/amu, an energy at which nuclear interactions leading to gamma ray 
lines also occur. The processes responsible for the X-ray transition are 
excitation by ambient gas particles and charge exchange to excited states of 
the iron ion. Then using nuclear reaction cross sections (for the most 

'recent values see Ramaty, Kozlovsky and Lingenfelter 1979), we calculated 

12 

the expected ratios of the production of 4.4 MeV line photons from C and 

56 

0 t 85 MeV photons from Fe to that of the 6.8 keV iron line. The results 
are shown in Figure II-6, 

We applied these results to possible observations of a 4.4 MeV line 
from the directions of the galactic center and the radio galaxy Centaurus A. 
For the former, if the line was produced by diffuse fluxes of low energy 
cosmic rays, the X-ray data can be made consistent with the gamma ray 
observation, but only if the energy in fast particles is so large that it 
conflicts with the vertical distribution of gas in the nuclear region. 

Possible resolutions of this conflict are obtained by confining the 
energetic particles to regions that are either hotter than ~ 10^K, or dense 
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enough to be opaque to the X-rays (by photoelectric absorption), Another 

possibility is time variability on the scale of years; this time scale would 

be consistent with the existence of a massive black hole at the galactic 

center. Such an object has been proposed to account for the time variability 

seen in Cen A. In addition, we find that the X— ray data from Cen A is 

inconsistent with the reported width of the 4,4 MeV feature there, assuming 

production by low energy cosmic rays. 

Also in Chapter II,. we have investigated the possibility of observing 

very narrow gamma ray emission lines from nuclei excited in interstellar 

dust grains. The densities of such grains are many orders of magnitude 

larger than the interstellar gas, so large that excited nuclei can be 

stopped before emitting the ga mma ray. Since the line widths are determined 

by Doppler broadening, a very narrow component of the line could result. 

Three conditions should be met: (i) a significant fraction of the element 

-must be condensed in grains, (ii) the lifetime of the nucleus* excited state 

—13 

must be long enough to allow deceleration Q> 10 s) , and (iii) the grains 
must b'e large enough that the recoiling nucleus does not pass out of it 
before being decelerated (linear dimensions ^ few tenths of a micron) . 

Figure II-8 through 11-11 show that this effect is important for several 
elements; thus observations of these very narrow components will yield 
information on the composition and size distributions of interstellar grains. 
We point out that this effect can only be observed with high resolution 
detectors. Future calculations should take better account of the grain 
composition and deviations from spherical symmetry. 

Spectral features in the hard X-ray region resulting from strong magnetic 
fields were considered in Chapter III. There we confirmed the results of 
Daugherty and Ventura (1977 and 1978) for the cyclotron line absorption 
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and emission rates, using a different representation of the electron wave 
functions,’ and we also calculated the relativistically correct electron- 
ion Coulomb cross sections. The momentum transfer in this latter process 
is quantized by the field, and significant modifications to the non-magnet ized 
case are apparent. Comparing the rates for these processes in the densities 
assumed to exist in the polar caps of the binary X-ray pulsars led us to 
conclude that electrons spend most of their time in the magnetic- ground state. 
In other words, cyclotron radiation is so efficient compared to collisional 
excitation that the perpendicular electron temperature is essentially zero. 

Energetic ions accreting onto an X-ray pulsar in a binary system are 
thought to be channeled into the polar cap with zero pitch angle. We have 
■used the above calculations to determine their energy loss rates and 
cyclotron line multiplicities, and the rate for production of cyclotron 
lines in a one- dimensional thermal gas. 

These results were applied to recent observations of a possible line 

feature in the spectrum of Hercules X— 1 at about 60 keV. Under either 

20 —3 

^seunjption, we found that electron densities are high enough (n >10 cm ) 
to provide large opacity in the line, as calculated from the absorption 
cross section. These considerations showed that the narrow reported 
width of the feature implies that the line emission is highly beamed at 
angles close to 90° with the polar cap field lines (a fan beam). 

Taken together with, the probable existence of an Alfven shell, our calcula- 
tions make the interpretation of a cyclotron emission line unlikely. The 
data are consistent with an absorption feature at lower energies, and we 
have presented a possible source of absorbing electrons, the knock-ons 
produced -T>y the infalling ions. Again, we turn to high resolution detectors 
for discriminating between the possibilities. In addition, further calcula- 

f * 

tions are needed to determine the electron and photon spectra in a self- 



consistent manner. We are presently considering the treatment of this 
» 

problem by means of a Monte Carlo simulation. 

Chapter XV treats in detail the annihilation of positrons in the 
galaxy and the resulting emission of 0.511 MeV gamma rays. By comparing 
their energy losses to the cross section for forming positronium by charge, 
exchange with neutral hydrogen, we have first evaluated the fraction- that 
annihilate before thermalizing. Since this cross section has a threshold 
at 6,8 eV, positrons which lose energy fast enough will fall below this 
threshold and will thermalize before annihilating. Energy losses are faster 
in an ionized gas, and we have found that the ionization fraction must be 
less than ~'5% for a majority of positrons to annihilate in flight. Then 
we have determined the rates for the various annihilation channels of 
thermal positrons. Positronium formation in the triplet state leads to 
three photon decay, while formation in the singlet state 'gives the two 
0.511, MeV photons upon annihilation, as is-the case for direct annihilation. 

Next we evaluate the shape of the radiation resulting from all, the 
processes considered in order to compare the' calculations with an observed 
line from the direction of the galactic center. The observed width requires 
a temperature of 10 5 K, and an ionization fraction > 5%. Because most of 
the matter in this region is in cold molecular clouds with practically no 
ionization, it is difficult to reconcile the observations with our calcu- 
lations unless the positron source is localized near the center or unless 
they are somehow' excluded from the clouds. We note that the first possibility 
is, consistent with the existence of the massive black hole mentioned in 
connection with the results of Chapter II. Among the various galactic 
regions, the warm component of the interstellar medium, filiments or knots 
in supernova remnants, and the' extended HII- region in the nuclear disk 
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could be sites for forming the 0.511 MeV line from positron annihilation, 

The calculated profiles show again that high resolution detectors can provide 
valuable information on the annihilation regions. Future calculations to be 
undertaken include following time variability of the sources and accounting 
for helium and heavier elements in the charge exchange and energy loss 
processes. 



APPENDIX • 


In this appendix, we present the details of the calculations of the 
cross sections and emission rates presented in Chapter III. First we eval- 
uate an integral occurring in' all the calculations: 

in an operator representation. The argument of the exponential can be 
written using equations III-8, 

lk x '*± * L k + (x~+ ) f.t k~(X*- ~ ) (.4 -a) 

where a vector $ is decomposed according to ■ 


(V*.±£ V») 


(/I -3) 


The commutation relations among the various operators are as follows: 
tt'3_ 

[)( + , X ]_ * ^ (/i -»0 

I**, ' TT± 'l- 0 

From Her zbacher, (1970, p. 167) we obtain for this situation: 

4*' + *' - £ Aj e A * e ~i 


Applying this relation to the exponential in A-l gives: 




krfr* - * 1 ?' 


■4- - £' ~ ' -t £ ^ C- 






(/?-0 


F Z 

where 5 = kj^ /(2eB). Thus, substituting in A-2 and using the hermicity of 
the operators, (JT*> = X. ; (tt ) - tr , yields 


f lkrx*-£*’ 


Expanding the exponential and operating on the functions according to (1II-9) 


k j.«-) (*v 


n 
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and (I 11-10) gives 

t - e' 3 

- L *.vvk- 


Ki* « 

J 1 *"/■>•• r ' ‘ (*-,*) 

which yields after using the orthogonality property of the f’s and 
k^ = kj_ e^^/y^ , where cp is the azimuth of It in the transverse plane, 

in,! „a J,, 

. • £ . C^S y^ n ' n^- V^ (ru-n, • C^) 

V*«*G»>vn,) U,-^)! iy^7)L 

Substituting T] = T^-mak (o, n^-n,^) and V- ~ (o, q.j-q 2 ) y iel(is 

" e ~ ? l ix ^ A ' 1M,u "^( ni ! pi 
^ t, 7Z fvo 


/* 

From the definition of the Laguerre polynominal 


(a- to) 


, t-0 . 

roTo ■ 05 ' (“ ml! P 

we see that (A-10) can be written 


J C-*) ™ p (.i-rc{4 | ) 
*— — o) 1 / 1 "’vTTTl ' 


(/Hi ) 


I n ,v^ct0 - ^ ««.(, )! 7 ' 

^ 3 ; C <3) £ f> 


^-li) 


We can now proceed to the evaluation of the cross sections. 

/ First, we calculate the absorption cross section and emission 

rate for line photons as prescribed by equations (III-13) and III-14) 

s fl ie (§?)*■ W* T* tp : <.*) (/^-J3') 

for absorption. Substitution from equation (XII-6) for tjr gives 

4,1 ' ^ Ww”w; uP5t^) ) ^ j' 4 "' 1 
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which, after using the hermiticity of the Dirac matrices (y^ 1 = y°y V ’y°)y 


becomes \ V f ». ■ . 

s f i * l~z 'r ' 1 

e -« <«»-*.*> u. 

This expression can be reduced to 2x2 from 4x4 by using 

I'D, > - U:) &-.») 

, / S, i \ -* 

where }. s y- = [ 1 “* ) an< ^ cr are the 'Pauli matrices. After commuting the expo - 

> 

nential factor to the left, we obtain the matrix operator between the states 


* c I** ( VV-Vv*,] V't Cv T ir y -»*0 (a - i 7 ) 

__ j tlTo + Ul + l^eX? * n +i £ *r)4 CtT° ti tr- fl TLa. k. K i 

\ /W ^ J 

where some components are not shown because they are not needed. Then 


the matrix element A-15 becomes 


£Oo-+a*»e-io3 £.n frL«r.£^n + £tn,<»n«) s? 't^'l| *:){■ h} 

J i / Jt 3 i ii Jl 

This equation is valid for all transitions between levels. We 
simplify by requiring that the initial state he the ground state, by letting 
q^ = o without loss of generality, and by evaluating in the center of 
parallel momentum frame, where » o and W^ = W^ s / a a +2eB. Under these 
conditions, equation A-18 becomes 

5. - U u/‘4a)^,)$ c^) 5 r r . 

fL aXt ItfVn/jWfc ViL U ^ + ^ +w )^^'C.w^mKK'f^A + r3j 

^ f 4 -i £“ . o fQ + £ L 'p. RWi+^XkV-fc*W) 

-^£' f; ]x w>o , 0 ctoj (/'-'‘rt 
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We note that if 0 is the polar angle of k and cp its azimuth, 

S ( .^ o) 

wher§ e e.^ are polarization- amplitudes in the ordinary and extra- 
ordinary modes. Substituting for the I's from A-12 gives 

f JtL £.«*» v * 5 3 '*■ ^ '“Tff 7* / ^r-1 JO^+av^ 


/t Ojvw,: *i, Luii+MAw, mu)}'*-^’ ' \P^T l^-L ^ J ' e3 ‘ 

*■ M«\coS& ffl + + ^ f ,t ^CCW^ « t ^U) c«*&) Z.O ft LvJc*<* e >u)) C0S& ix ] j 1 (fit -Zl) 

Carrying out the same analysis for the emission case gives (assuming p., q. = 0) 

r T't & (~tfSe iV J %f f r ,f ^ 

•Jft. — ~ 2 L r77Tr777 _ TT7 7 r . . t £ 7 *•_ J \ f CtT* 1 Ifm 


f': - £& / •; r r# 

/i. - <*l* Wj W f (vO t vm t H^ + »k)] Vl ‘ € ffy? I 

f ^ f^)] f u) £( W^tn^+WC^Ofo-t (.Yif 4 W e +U)) <^S$- J J (/Wl) 

The absorption cross section is obtained by squaring the magnitude of 
and summing over electron final states. The square of the 6 -functqons are 


handled heuristically in Bjorken and Drell (1964): 

1 *C^| = \Scw)p= S 


0-0 


ttfxere L and T are the (large) normalization length parallel to the field 


and the normalization time. Then the rate of absorptions is 


7 <C>ta = y 


£ K-r 


(ft’**) 


where V is the electron normalization volume, <7 is the absorption cross 
section, and v ^ is the relative velocity of the photon and electron. 
Solving for the absorption cross section yields equ,t±on (III-15) in the 
center of parallel momentum frame. 

The emission rate is shown in here in more detail to show the 
processes involved in arriving at equation III-15, First, the sum over 


q^ yields 


(±Use c r) 

~ffT 


c^if [* 


y* w 


tl f - ^-3 
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The Integral over is ignored because of the 6 function. Then the emission 
3 

rate into d k centered on k for electrons deexciting from the first excited 
state can be written as 

0«) 




(an)' T 

3 3 

where Vd k/ (2 tt) expresses the proportionality of the rate to the phase 
Space volume for the emitted photon. Carrying out the operations yields 

= w uTvffc/j f C^+m r +i0)5 + 

+ ^ri ^F fh 1 e) + tc^ + »««£.-+ U) CO S l fi) f 0 l ’+C\jy 4 + m < -tio) i C^ l & ^ 




When carrying out the integral over d k to obtain the inverse of the life- 
time T^» it must be noted that depends on to so that the § function in 
A-27 should be written 


Then we obtain 


-I _ — 


-f w- = ,- 7 — fc — S l, A- "^Wfr+m,W/P \ 




i 


jje 




j fz^jZi- 

where the braces represent the quantity in braces in A-27 . This integral 
was evaluated numerically; the results are shown in Figuee III-l. 

Next we consider the scattering of an electron by an ion in 
an external field, first order approximation, the relativistic extension of 
the Bom approximation. The matrix element involved in the ion’s rest frame 
is given by equation 111-20; we Fourier expand the ion potential 
Ze/r‘e 1 Sc , where r sc is the screening length: 


d» = /(0 e a,i 4, <ib , 


/U*J» 


<nr£ c. 

V * r£~ 


(4-30) 


Then equation III-20 becomes 

S ft ~ “ ,1rt J (tVr,; 1 ) 1 j<iV ■fy ^ tf* Lx) 


(A-*0 
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Substituting for 1} from equation III-6 and using hermiticity gives 

|r^ U<* -f * W 

v J %f 

Us + * L kf* «V +m *) 7 v e - * (>% niij (/j-5>2.) 


4 -4 

ik-x 


After moving e to the left,’ we reduce the inner 4x4 to the 2x2 repre- 
sentation as ’.in. A-16 through A-18: 


u 5 + i c [r° rr 0 - ?- uf4-£)+m e ] *° *rO u } - 

— ^ S f j [ CT^tyn ^ 1 4 - tt 4 k.»n + i cr- ltt+ Vl) *77 { r i 


(V**) 


^ 2_ •-) 0 — ^ 

We note that n x n = j V x (eA) = — B. Then substituting this expression 

into A-32 and assuming the initial state is the ground state gives: 
c ... .-,-i 


^ ^ + **0 H+ «c)] x J ( k\r 3 i x )' 1 ^ . t «ct> £ ‘ ^ ^ 

The elasticity of the interaction is expressed by the fact that 
the Coulomb potential had no time varying component. This results in a 
factor 5(W^-W^) in the above time integration. In evaluating the integral, 
we put = W to obtain: 

s K - Us- \ 0 . sUH<y--) -I,,- to 

{Vi [^WCW-fnuV^ - Vi (^-3r) 

Substituting for the I function and carrying out the k„ and cp integrations 


_2 

gives (ignoring r gc ): 


5* = Trijjc*- (“ iS ' Ivvl f £ 

^ L Vj CW + nO Jo 5 *■ ((>t'K)Vc3^ 8 ^ 'J* 

4 if Vi CWjkm ^ ~ V». ^ • 

We rewrite the above, noting = q +j f and putting j = j f 

*- v mi Uii&n^i******* 



y^ip 


(^- 37 ) 



7.9 


-TV' 

£; £S 
* & 


where x = (p^-p^) /(2eB). Next we square this expression, sum over final q 

and average over initial q according to l/Q § where Q is obtained by 

*ii=o 

noting that the distance from the guiding center to the field line q is 


given by 


tt’lv (xVfx-xOJ - ft- ViH* . 

2 

The maximum q, Q, in a given area A(= ttR ) is given by 


<X - Q+i. " i ei3R T 


The rate- is obtained from 


- cl 


(W-ss - ) 


(A- 3<0 


*4 W - i <f W lvt N ,| - t 1 I I VI 1 -Vo) 

In the square, the 8 function is handled as in A-23, The 8-function remains, 
and is used in the integral, after 

solving (A-40) for. cr^ (noting that v = p/W) gives 

j ^ ±" l»°-M f\-£ i 

00 p « k> • - 12 . 0-1 ' c) J 

\ ^ f 1 . m^r~ J (A -41) 

The sum can be rewritten by using first of all the Laguarre expansion theorem 


(Danese, 1965, p. 517): 

61.%') - Z- o Sj" 1 Li 

Allowing G(S) = (5+x) and a = j • gives the expression 


(A-«> 


e. rr«£ <**»>)«$ f 


But by A-42, we have 


1 ?->* if- i 1 u v s >- h 


0 )-^) 


Finally we obtain 

cr( P -)= JllJk 2 ' _ K ^ 


I^Kvi^y & rif 
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2 2 

where r = e /(.me ) is the classical electron radius, and the e are 
e J 

given by (putting 5+x‘ = t in A-45) 
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TABLE II-l 


Final 


Element 

Abundance 

State 

E = 5 MeV/amu 

10 

13 

20 

30 

H 

1 . 

Is 

2.73-21 

8.16-23 

2.24-23 

2.41-24 

3.08-25 



2p 

2.25-21 

6,64-23 

2.00-23 

1.55-24 

1.58-25 



2s 

7,92-23 

3.58-24 

1,00-24 

1.12-25 

2.25-26 



Is 

4.94-21 

1.53-22 

4.21-23 

4.57-24 

5.89-25 

He 

0.085 

2p 

4.08-21 

1.24-22 

3.11-23 

2.94-24 

3.02-25 


(0.1) 

2s 

1.46-22 

6,72-24 

1.90-24 

2.14-25 

4.29-26 


3.71-4 

Is 

4.54-21 

2.15-22 

6.48-23 

8.14-24 

1.15-24 

C 

( 2 , 1 - 3 ) 

2p 

3.68-21 

1.72-22 

4.79-23 

5.11-24 

5.73-25 



2s 

1.78-22 

1.00-23 

3.05-24 

3.91-25 

8.63-26 


1,18-4 

Is 

1.81-21 

1.01-22 

3.12-23 

4.11-24 

5.97-25 

N 

(1,3-3) 

2P 

1.48-21 

7.89-23 

2.27-23 

2.56-24 

4,57-26 



2s 

7.86-23 

4.72-24 

1.49-24 

2.01-25 

4,57-26 


6.76-4 

Is 

1.91-20 

1.22-21 

3.93-22 

3.52-23 

8.30-24 

0 

(3.5-3) 

2p 

1.56-20 

9.59-22 

2.86-22 

3.40-23 

4.08-24 



2s 

9.15-22 

5.84-23 

1.91-23 

2.74-24 

6.46-25 


1.08-4 

Is 

3.59-21 

3.00-22 

1.05-22 

1.67-23 

2.73-24 

Ne 

(1.9-4) 

2p 

2.90-21 

2.30-22 

7.47-23 

1.00-23 

1.29-24 



2s 

2.11-22 

1.45-23 

5.19-24 

8.64-25 

2.20-25 


3.34-5 

Is 

1.46-21 

1.50-22 

5.63-23 

1.01-23 

1.82-24 

Mg 

(2.0-4) 

2p 

1.17-21 

1.13-22 

3.93-23 

5.92-24 

8.79-25 



2s 

9.67-23 

7.03-24 

2.80-24 

5.59-25 

1.54-25 


3.14-5 

Is 

2.25-21 

2.55-22 

1.03-22 

2.07-23 

4 .09-24 

Si . 

(2.3-4) 

2p 

1.77-22 

1.80-22 

6.96-23 

1.16-23 

1.80-24 



2s 

1.40-22 

1.12-23 

5 . 05-24 

1.23-24 

3.61-25 


1.57-5 

Is 

1.35-21 

1.51-22 

6.42-23 

1.42-23 

3.08-24 

S 

(8,0-5) 

2p 

1.06-21 

1.07-22 

4.12-23 

7.60-24 

1.28-24 



2s 

6.76-23 

7.03-24 

3.38-24 

9.15-25 

2.83-25 


2.61-5 

Is 

6.75-21 

6,62-22 

2.85-22 

6.51-23 

1.75-23 

Fe 

'(1,6-4) 

2P 

4,72-21 

4.55-22 

1.61-22 

2.82-23 

5.52-24 



2s 

3.89-22 

3.59-23 

1.74-23 

4.78-24 

1.72-24 



Is 

4.85-20 

3.29-21 

1,19-21 

2.01-22 

4.02-23 

TOTALS 


2p 

3.87-20 

2.49-21 

7.94-22 

1.09-22 

1.62-23 



2s 

2.30-21 

1.59-22 

6.22-23 

1.20-23 

3.58-24 



TABLE 11-2 


Multiplicities 


E(MeV/amu) 

Neutral Medium 
(Photons) 

Ionized Medium 
(Photons) 

2 

-4 • 

2.0x10 

2,4xl0~ 5 

4 

.0296 

,00414 

6 

.281 

.0413 

8 

1.031 

,156 

10 

2.32 

.361 

12 

3.92 

.620 

14 

5.48 

.877 

16 

6,77 

1.09 

18 

7.75 

1.26 

20 

8.46 

1,38 

22 

8.98 

1.47 

24 

9.38 

1.54 

26 

9,68 

1.60 

28 

9.92 

1.64 

30 

10,11 

1.67 

32 

10.26 

1,70 

34 

10.38 

1.72 

36 

10.48 

1.74 

38 

10.56 

1,76 


40 


10.63 


1.77 





PArtE 


Photon 

Energy 

(KeV) 

Emission 

Mechanism 

Production 

Processes 

'feaa life 
(sec) 

* Photon 
Energy 
(MeV) 

Emission 

Mechanism 

Production 

Processes 

Mean Life 
(see) 

0.092 

55^*1.608,55^*1 316 

' S W P a> 5S l*‘ 

- 

1.334 

52 Cr *2.763^52 Cr *l 434 

56 , .52 * 

Fe(p.x) Cr 

__ 

c* no 



8-SxlB" 10 

1.369 

2 ‘h R * 1 - 3 ®%.6. 

24 , ,"24„ * 

Mg(p.p / ) Mg 

1.75xl0' 12 

■ 1SS 

54 c*»- IS %. s 

5 ®Fe<p,n) 5 ®Go* 

- io-‘° 



24„ , ,.24., * 

Hgfcr.a ) Mg 

].7Sx!ir 12 

0 197 

19 r*° 19 s 

2 V(p.2?) ,9 F* 

1 .3xl0 -7 



•>8 74 * 

“siCp.x)^ Kg 

1.75xl0' 12 

0.238 

, V‘°- 23 Vx. 

2 °He(p.pn) 19 Pe 

2 6xl0 -8 

1.370 

“f.* 2 - 381 . 5 ^* 0 931 

S ®Fe(p,pn) SS Fe* 



0.275 

19., *0.275 
Ne -»g.s. 

20 Ne(p,pn) l9 Ko* 

6.lxl(f n 

1.408 

"fp* 2 - 408 ^.*. 

S6 Fie(p.x) 54 Fc* 

l.JxIO* 12 

0.412 
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